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Abstract 

Recent experimental and theoretical work has shown that there are condi- 
tions in which a trapped, low-density Bose gas behaves like the one-dimensional 
delta-function Bose gas solved years ago by Lieb and Liniger. This is an intrin- 
sically quantum-mechanical phenomenon because it is not necessary to have a 
trap width that is the size of an atom - as might have been supposed - but it 
suffices merely to have a trap width such that the energy gap for motion in the 
transverse direction is large compared to the energy associated with the motion 
along the trap. Up to now the theoretical arguments have been based on vari- 
ational - perturbative ideas or numerical investigations. In contrast, this paper 
gives a rigorous proof of the one-dimensional behavior as far as the ground state 
energy and particle density are concerned. There are four parameters involved: 
the particle number, N, transverse and longitudinal dimensions of the trap, r 
and L, and the scattering length a of the interaction potential. Our main re- 
sult is that if r/L — > and N — > oo the ground state energy and density can 
be obtained by minimizing a one-dimensional density functional involving the 
Lieb-Liniger energy density with coupling constant ~ a/r 2 . 

This density functional simplifies in various limiting cases and we identify 
five asymptotic parameter regions altogether. Three of these, corresponding to 
the weak coupling regime, can also be obtained as limits of a three-dimensional 
Gross-Pitaevskii theory. We also show that Bose-Einstein condensation in the 
ground state persists in a part of this regime. In the strong coupling regime the 
longitudinal motion of the particles is strongly correlated. The Gross-Pitaevskii 
description is not valid in this regime and new mathematical methods come into 
play. 

*Work partially supported by U.S. National Science Foundation grant PHY 01-39984. 
^Erwin Schrodinger Fellow, supported by the Austrian Science Fund. 
© 2003 by the authors. This paper may be reproduced, in its entirety, for non-commercial purposes. 
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1 Introduction 

The technique of trapping and cooling atoms, that led to the first realization of 
Bose-Einstein condensation (BEC) in dilute alkali gases in 1995 [5, 18], has recently 
opened the possibility for experimental studies, in highly elongated traps, of Bose 
gases that are effectively one-dimensional. Some of the remarkable properties of 
ultracold one-dimensional Bose systems with delta function interactions, analyzed 
long ago [22, 23], may thus become accessible to experimental scrutiny in the not too 
distant future. Among these are pseudo-fermionic behavior [12], the absence of BEC 
in a dilute limit [21, 36, 14, 34, 10], and an excitation spectrum different from that 
predicted by Bogoliubov's theory [23, 17, 20]. The paper [33] by Olshanii triggered 
a number of theoretical investigations on the transitions from 3D to an effective 
ID behavior with its peculiar properties, see, e.g., [6, 7, 8, 11, 13, 19, 32, 35, 39]; 
systems showing the first evidence of such a transition have recently been prepared 
experimentally [4, 15, 16, 38]. 

Until now the theoretical work on the dimensional cross-over in elongated traps 
has either been based on variational calculations, starting from a three-dimensional 
delta-potential [7, 13, 33], or on numerical Monte Carlo studies [1, 3] with more re- 
alistic, genuine 3D potentials but particle numbers limited to the order of 100. This 
work is important and has led to valuable insights, in particular about different pa- 
rameter regions [8, 35, 32], but a more thorough theoretical understanding is clearly 
desirable since this is not a simple problem. In fact, it is evident that for a potential 
with a hard core the true 3D wave function does not approximately factorize in the 
longitudinal and transverse variables and the effective one-dimensional potential can 
not be obtained by simply integrating out the transverse variables of the 3D poten- 
tial. In this sense the problem is more complicated than in a somewhat analogous 
situation of atoms in extremely strong magnetic fields [2, 30], where the Coulomb 
interaction behaves like an effective one-dimensional delta potential when the mag- 
netic field shrinks the cyclotron radius of the electrons to zero. In that case the delta 
potential can be obtained formally by integrating out the variables transverse to the 
field in a suitable scaled Coulomb potential. With a hard core, on the other hand, 
where the energy is essentially kinetic, this method will not work since it would 
immediately introduce impenetrable barriers in ID. The one-dimensional effective 
interaction emerges only if the kinetic part of the Hamiltonian and the potential are 
considered together. 

In the present paper we start with an arbitrary, repulsive 3D pair potential of 
finite range and prove rigorously that in a well defined limit the ground state energy 
and particle density of the system are described exactly by a one-dimensional model 
with delta-function interaction. This is a highly quantum-mechanical phenomenon 
with no classical counterpart, since a ID description is possible even though the 
transverse trap dimension is much larger than the range of the atomic forces. It 
suffices that the energy gap associated with the transverse confinement is much 
larger than the internal energy per particle. 

While the three-dimensional density remains low (in the sense that distance be- 
tween particles is large compared to the three-dimensional scattering length) the 
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one-dimensional density can either be high or low. We remark that, in contrast to 
three-dimensional gases, high density in one dimension corresponds to weak inter- 
actions and vice versa [22]. In this paper we shall always be concerned with large 
particle number, N, which is appropriate for the consideration of actual experiments. 
In order to make precise statements we shall typically take the limit N — ► oo but the 
reader can confidently apply these limiting statement to finite numbers like N = 100. 

Besides N, the parameters of the problem are the scattering length, a, of the 
two-body interaction potential, and two lengths, r and L, describing the transverse 
and the longitudinal extension of the trap potential, respectively. To keep the intro- 
ductory discussion simple let us first think of the case that the particles are confined 
in a box with dimensions r and L. The three-dimensional particle density is then 
p 3D = N/(r 2 L) and the one-dimensional density p 1D = N/L. The case of quadratic 
or more general trapping potentials will be considered later. We begin by describing 
the division of the space of parameters into two basic regions. This decomposition 
will eventually be refined into five regions, but for the moment let us concentrate on 
the basic dichotomy. 

In earlier work [27, 26] we have proved that the three-dimensional Gross- 
Pitaevskii formula for the energy (including its limiting 'Thomas-Fermi' case) is 
correct to leading order in situations in which a is small and N is large. This en- 
ergy has two parts: The energy necessary to confine the particles in the trap, which 
is roughly (h 2 /2m)N(r~ 2 + L~ 2 ), plus the internal energy of interaction, which is 
(h 2 /2m)N4-jrap 3I) . The trouble is that while this formula is correct for a fixed con- 
fining potential in the limit N — > oo with a 3 /? 3D — ► 0, it does not hold uniformly if 
r/L gets small as N gets large. In other words, new physics can come into play as 
r/L — > and it turns out that this depends on the ratio of a/r 2 to p 1D = N/L . As 
we shall show, the two basic regimes to consider in highly elongated traps, i.e., when 
r«L, are 

• The one-dimensional limit of the three-dimensional Gross-Pitaevskii/'Thomas- 
Fermi' regime 

• The 'true' one-dimensional regime. 

The former is characterized by aL/r 2 N — > 0, while in the latter regime aL/r 2 N is of 
the order one or even tends to infinity (which is referred to as the Girardeau- Tonks ^ 
region). These two regimes correspond to high one-dimensional density (weak in- 
teraction) and low one-dimensional density (strong interaction), respectively. The 
significance of the combination aL/r 2 N can be understood by noting that it is 
the ratio of the 3D energy per particle, ~ ap 3T> ~ Na/r 2 L, to the ID energy 
~ (p 1D ) 2 = {N/L) 2 . Physically, the main difference between the two regimes is 

^We call this the Girardeau- Tonks region only because many authors refer in the present context 
to Tonks [41]. In our opinion this should really be called the Girardeau region because it was he 
who first understood how to compute the spectrum of a ID quantum-mechanical hard core gas and 
who understood that the Fermi-Dirac wave functions played a role [12]. Tonks was interested in 
the positive temperature partition functions of a hard core classical gas - a very different and much 
simpler question. 
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that for strong interactions the motion of the particles in the longitudinal direction 
is highly correlated, while in the weak interaction regime it is not. Mathematically, 
this distinction also shows up in our proofs. 

In both regimes the internal energy of the gas is small compared to the energy 
of confinement which is of order N/r 2 . However, this in itself does not imply a 
specifically one-dimensional behavior. (If a is sufficiently small it is satisfied in a 
trap of any shape.) One-dimensional behavior, when it occurs, manifests itself by 
the fact that the transverse motion of the atoms is uncorrelated while the longitudinal 
motion is correlated (very roughly speaking) in the same way as pearls on a necklace. 
Thus, the true criterion for ID behavior is that aL/r 2 N is of the order unity and not 
merely the condition that the energy of confinement dominates the internal energy. 

The starting point for our investigations is the Hamiltonian for TV spinless Bosons 
in a confining 3D trap potential and with a short range, repulsive pair interaction. 
We find it convenient to write the Hamiltonian in the following way (in appropriate 
units): 

N 

H N ^ a = Y,(-^ + V r ± (xf) + V L (z 3 ))+ "-(l^-^-l) ( L1 ) 

j=l l<i<j<N 

with x = (x, y, z) = (x -1 , z), 

V r H^) = 4^(x», V L (z) = ±V(z/L), „ a (|x|) = 4^(|x|/a) , (1.2) 

where r, L, a are variable scaling parameters while V- 1 , V and v are fixed. The inter- 
action potential v is supposed to be nonnegative, of finite range and have scattering 
length 1; the scaled potential v a then has scattering length a. The external trap 
potentials V and V 1 - confine the motion in the longitudinal (z) and the transversal 
(x- 1 ) directions, respectively, and are assumed to be locally bounded and tend to oo 
as \z\ and jx- 1 ! tend to oo. To simplify the discussion we find it also convenient to 
assume throughout that V is homogeneous of some order s > 0, namely V(z) = \z\ s , 
but weaker assumptions, e.g. asymptotic homogeneity [28], would in fact suffice. 
The case of a simple box with hard walls is realized by taking s = oo, while the usual 
harmonic approximation is s = 2. Moreover, to avoid unnecessary technicalities we 
shall assume that V 1 - is polynomially bounded at infinity, but our results certainly 
also hold for faster growing potentials, or even finite domains with Dirichlet bound- 
ary conditions. Units are chosen so that h = 1 and 2m = 1. It is understood that 
the lengths associated with the ground states of —d 2 /dz 2 + V(z) and —A 1 - + y- L (x- L ) 
are both of the order 1 so that L and r measure, respectively, the longitudinal and 
the transverse extensions of the trap. We denote the ground state energy of (1.1) by 
.E;Q M (jV, L,r,a) and the ground state particle density by r a (x). 

In parallel with the three-dimensional Hamiltonian we consider the Hamiltonian 
for n Bosons in one dimension with delta interaction and coupling constant g > 0, 
i.e., 

n 

H ™ = Y.- d l + 9 E 8{zi- Zj ), (1.3) 

j=l l<i<j<n 
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where dj = d/dzj. We consider this Hamiltonian for the Zj in an interval of length 
£ in the thermodynamic limit, I — > oo, n — > oo with p = njt fixed. The ground 
state energy per particle in this limit is independent of boundary conditions and can, 
according to [22], be written as 



el D (p)= P M9/p) , 



(1.4) 



with a function e(t) determined by a certain integral equation. Its asymptotic form 
is e(t) « \t for t < 1 and e(t) -> 7r 2 /3 for i -> oo. Thus 



e 0°(p) ~ for 5/p < 1 



and 



7T 



eo (P) ~ yP ^ S/P» 1 • 



(1.5) 



(1.6) 



Taking pe^(p) as a local energy density for an inhomogeneous one-dimensional sys- 
tem we can form the energy functional 



£[P] = / (IVVP(^)I 2 + V L {z)p{z) + p{zfe{g/p{z))) dz 



(1.7) 



with ground state energy defined by minimizing £[p] over all normalized densities p, 
i.e., 



E(N, L, g) = inf \s [p] : p > , jf p(z)dz = ivj . 



(1.8) 



By standard methods (cf., e.g., [27]) one can show that there is a unique minimizer, 
i.e., a density pN,L,g{z) with J p^ ! L,g(z)dz = N and £ [pN,L,g] = E(N, L, g). Here it 
is important to note that t <— > t 3 e(l/t) is convex. We define the mean ID density of 
this minimizer to be 

P=Tj [ (PN,L, 9 {z)) 2 ^ . (1.9) 

In a rigid box, i.e., for s = oo, p is simply AT/L, but in more general traps it depends 
also on g besides iV and L. The order of magnitude of p in the various parameter 
regions will be described in the next section. 

Our main result relates the 3D ground state energy of (1.1), E® M (N, L,r,a), to 
the ID density functional energy E(N, L, g) for a suitable g in the large N limit 
provided r/L and a/r are sufficiently small. To state this precisely, let e 1 - and 
^(x- 1 -) denote the ground state energy and the normalized, nonnegative ground state 
wave function of — + V r_L (x^), respectively. The corresponding quantities for 
— A -1 + Vf L (x. ± ) are e^/r 2 and 6 r (x J -) = (l/r)6(x ± /r). In the case that the trap is a 
cylinder with hard walls b is a Bessel function; for a quadratic V 1 - it is a Gaussian. 
In any case, 6 is a bounded function and, in particular, b G L 4 (IR 2 ). Hence we can 
define g by 



87TO 
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I 



|&(x x )|Vx x = 8vra / |6 r (x x )|Vx 



/ 



(1.10) 



Our main Theorem is: 
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Theorem 1.1 (Prom 3D to ID). Let N — > oo and simultaneously r/L — > and 
a/r — > in suc/i a way i/ia£ r 2 p • min{p, g] — ► 0. T/ien 

E(N,L,g) K > 

Note that because of (1.5) and (1.6) the condition r 2 p- min{p, g} — > is the same 

as 

eJ D (p) « 1/r 2 , (1.12) 

i.e., the average energy per particle associated with the longitudinal motion should 
be much smaller than the energy gap between the ground and first excited state of 
the confining Hamiltonian in the transverse directions. (The precise meaning of <C is 
that the ratio of the left side to the right side tends to zero in the limit considered.) 
Note also that while the one-dimensional density can be either high or low (compared 
to g), the gas is always dilute in a three-dimensional sense in the limit considered, 
i.e., a 3 p m ~ a 2 gp < 1. 

The two regimes mentioned previously correspond to specific restrictions on the 
size of the ratio g/p as N — > oo, namely j/p < 1 for the limit of the 3D Gross- 
Pitaevskii regime (weak interaction/high ID density), and g/p > for the 'true' 
one-dimensional regime (strong interaction/low ID density). We shall now describe 
briefly the finer division of these regimes into five regions altogether. Three of them 
(Regions 1-3) belong to the weak interaction regime and two (Regions 4-5) to the 
strong interaction regime. In each of these regions the general functional (1.7) can be 
replaced by a different, simpler functional, and the energy E(N, L,g) in Theorem 1.1 
by the ground state energy of that functional. 

The five regions are 

• Region 1, the Ideal Gas case: g/p <C iV -2 , with p ~ N/L, corresponding 
to a non-interacting gas in an external potential. 

• Region 2, the ID GP case: g/p ~ N~~ 2 , with p ~ N/L, described by a ID 
Gross-Pitaevskii energy functional with energy density \gp 2 ■ 

• Region 3, the ID TF case: iV~ 2 < g/p <C 1, with p ~ (N/L)(NgL)- 1 ^ s+1 \ 
where s is the degree of homogeneity of the longitudinal confining potential V. 
This region is described by a Thomas-Fermi type functional with energy density 
\gp 2 ■, without a gradient term. 

• Region 4, the LL case: g/p ~ 1 , with p ~ (N/ L)N~ 2 /( s+2 \ described by an 
energy functional with the Lieb-Liniger energy (1.4), without a gradient term. 

• Region 5, the GT case: g/p > 1, with p ~ (N / L)N- 2 ^ s+2 \ described by 
a functional with energy density ~ p 3 , corresponding to the Girardeau- Tonks 
limit of the LL energy density. 

We note that the condition g/p ~ 1 means that Region 4 requires the gas cloud to 
have aspect ratio r/L of the order N^ 1 (a/r) or smaller, where L = N/p ~ LN 2 ^ S+2 ^ 
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is the length of the cloud. Experimentally, such small aspect ratios are quite a 
challenge and the situations described in [4, 15, 16, 38] are still rather far from this 
regime. It may not be completely out of reach, however. 

The condition a/r — > is automatically fulfilled in Regions 1-4, provided (1.12) 
holds, since (a/r) 2 = r 2 (gp)(g/p) = r 2 p 2 (g/p) 2 <C 1 by (1.12) if g/p is bounded. 
Moreover, as discusssed in the next section, in Regions 1-2 the condition r/L — > 
implies (1.12) and hence a/r — > , and in Region 4, a/r — > implies (1.12). The 
hypotheses of Theorem 1.1 are thus not entirely independent. 

In the next Section 2 we define the various energy functional more precisely and 
also discuss the ID behavior of the density PjVLra( x )' separately for each region. 
Moreover, we prove, in Subsection 2.7, that Regions 1-3 can be reached as limiting 
cases of a 3D Gross-Pitaevskii theory. In this sense, the behavior in these regions 
contains remnants of the 3D theory, which also shows up in the fact that BEC prevails 
in Regions 1 and 2, as discussed in Section 5. Heuristically, these traces of 3D can 
be understood from the fact that in Regions 1-3 the ID formula for the energy per 
particle gp ~ aN/(r 2 L), gives the same result as the three-dimensional formula [31], 
i.e., scattering length times three-dimensional density. This is no longer so in Regions 
4 and 5 and different mathematical methods are required. 

Despite significant differences the proof of Theorem 1.1 has some basic strategies 
in common with the (considerably simpler) proof of the Gross-Pitaevskii limit Theo- 
rem in [27] (see also [31] and [26]). The upper bound for the energy in Regions 1-3, 
given in Subsection 4.1, is the simplest estimate and can be obtained by a method 
analogous to that of [27]. For the lower bound, and also for the upper bound in Re- 
gions 4-5, one considers first finite numbers of particles in boxes with Neumann or 
Dirichlet boundary conditions, and subsequently puts these boxes together to treat 
inhomogeneous external potentials and the infinite particle number limit. For the 
lower bound of the energy in the boxes Dyson's Lemma [9, 31, 26], which converts a 
"hard" potential into a "soft" potential at the expense of sacrificing kinetic energy, is 
an essential tool. The main differences compared to [27] are on the one hand due to 
the fact that in Regions 1-3 the lower bounds in [27] are not valid because they are 
not uniform in the shape of the trap, and on the other hand due to the correlations 
in the Lieb-Liniger wave function for the longitudinal motion in Regions 4-5. For 
the latter reason even the proof of the upper bound in Regions 4-5 is considerably 
more involved than for Regions 1-3. 

Section 3 contains the main technical estimates for the boxes. We consider here 
the 3D Hamiltonian for a finite number of particles on a finite interval in the lon- 
gitudinal direction with Neumann or Dirichlet boundary conditions and a confining 
potential in the transverse directions. We estimate its energy from above and be- 
low in terms of the energy of a ID Hamiltonian with delta interactions. In Section 4 
we apply these results to prove Theorem 1.1. In the last Section 5 we consider the 
question of Bose-Einstein condensation and prove that it holds in Regions 1 and 2. 

Part of the results of this work were announced in [29] . 
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2 The five parameter regions 

We shall now discuss the simplifications of Theorem 1.1 in the five different parameter 
regions. Besides the ground state energy we also consider the convergence of the 
quantum mechanical density PyvLra( x ) avera g e d over the transverse variables, i.e., 

We state the results in the form of five theorems, one for each region. These theorems 
will be proved in Section 4. Together they imply Theorem 1.1 because the energy 
functionals involved are all limiting cases of the general functional (1.8). The methods 
for proving the latter statement are fairly standard (cf., e.g., [27, 28] for similar 
computations). Since a proof of all five limit theorems for the functional (1.8) would 
be largely repetitious, we shall limit ourselves to giving a proof of two of them in 
Subsection 2.6, as an example. 

2.1 The Ideal Gas Region 

This region corresponds to the trivial case where the interaction is so weak that it 
effectively vanishes in the large N limit and everything collapses to the ground state 
of —d 2 /dz 2 + Vl(z). By scaling the ground state energy and wave function of this 
latter operator can be written as L~ 2 e" and L~ l l 2 yf pW (z/L) where e" and ypHz) 
are the corresponding quantities for —d 2 /dz 2 + V(z). 

Theorem 2.1 (Ideal gas limit). Suppose r/L — > and NgL ~ NaL/r 2 — > as 
N -> oo. Then 

{N/L 2 )- 1 (e qm {N, L, a, r) - Ne ± /r 2 ^j -► e 11 (2.2) 

and 

(N/Ly'p^ljL^^pHz) (2.3) 

weakly in L 1 (IR) . 

2.2 The ID Gross-Pitaevskii Region 

This region is described by the ID GP density functional 

£l>] = / (|Vv^(z)| 2 + V L {z)p{z) + \gp{zf) dz (2.4) 

JR 

corresponding to the high density approximation (1.5) of the interaction energy in 
(1.7). Its ground state energy 

E GP (N,L,g) =inf |fgj[p] : p > , jf p{z)dz = ivj (2.5) 

has the scaling property 

E GP (N, L, g) = (N/L 2 )E GP (l, 1, NgL) (2.6) 
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(2.7) 



Theorem 2.2 (ID GP limit). Suppose r/L — > and NgL ~ NaL/r 2 is fixed as 
N -> oo. T/ien 



weakly in L 1 (IR) . 

Remark. If r/L — > and iVgL stays bounded, as in Regions 1 and 2, condition (1.12) 
is automatically satisfied, because gpr 2 ~ aN/L ~ (r / L) 2 (NgL). Likewise, a/r — > 0, 
because a/r = (r / L)N~ 1 (N gL). 

2.3 The ID 'Thomas-Fermi' Region 

This region is a limiting case of the previous one in the sense that NgL ~ NaL/r 2 — > 
oo, but a/r is sufficiently small so that g/p ~ (aL/Nr 2 )(NaL/r 2 ) 1 ^ s+1 ^ — > 0, i.e., 
the high density approximation in (1.5) is still valid. Here s is the degree of homo- 
geneity of V and the explanation of the factor (NaL/r 2 ) l /^ s+l ^ ~ (iVpL) 1 /^ 4 " 1 ) is as 
follows: The linear extension L of the minimizing density p^i g is for large values of 
NgL determined by Vl(L) ~ g(N/L), which gives L ~ {NgL) l ^ s+1 ^L. In addition 
condition (1.12) requires <C r -2 , which means that (Na/L)(NaL/r 2 y 1 ^ s+1 ^ -► 0. 

If iVgL ~ NaL/r 2 — > oo the gradient term in the functional (2.5) becomes 
negligible compared to the other terms. In fact, by a simple scaling, 




(2.8) 



and 



(2.9) 



£2 P 9 [P] = ^(NgLyf^ ^{(NgL) 



(2.10) 



where the scaled density p is determined by 



p(2) = (N / Ltf) p{z / Ltf) , with L TF := (iVgL) 1/(s+1) L 



(2.11) 



This leads to the functional 




(2.12) 



whose ground state energy 




(2.13) 



has the scaling property 



E TF (N,L,g) = (N/L 2 )(NgL) s ^ s+1 ^E TF (1,1,1) . 



(2.14) 
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The minimizer pJf L satisfies 

P T n F l, 9 (L tf z) = (N/L^pJl^z) (2.15) 
and can be computed explicitly: 

pJJ,i(z) = [p TF -V(z)} + (2.16) 

where [t]+ = max{i, 0} and // TF is determined by the normalization J pJ F 1 (z)dz = 1. 

Because of a formal similarity with the Thomas-Fermi energy functional for 
fermions, which also has no gradient terms, the functional (2.12) is in the litera- 
ture usually referred to as a 'Thomas-Fermi' functional, but the physics is, of course, 
quite different. 

The limit theorem in this region is 

Theorem 2.3 (ID 'TF limit'). Suppose, as N -> oo, r/L -> 0, NgL ~ NaL/r 2 -> 
oo, but g/p ~ (aL^p/Nr 2 ) — > and g/?r 2 ~ Na/L^p — > wii/t Ltf given in (2.11). 
Then 

{N/L 2 )- 1 {NgL)- s / { - s ^ |"^Q M (iV, L, a, r) - iV e ^/ r 2 j £ TF (1, 1, 1) (2.17) 

and 

{N/L^pXl^L^z) - ^(z) (2.18) 

weakly in L 1 ) 



2.4 The Lieb-Liniger Region 

This region corresponds to the case g/p ~ 1, so that neither the high density (1.5) 
nor the low density approximation (1.6) is valid and the full LL energy (1.4) has 
to be used. The extension L of the system is now determined by Vl(L) ~ (N/L) 2 
which leads to Z LL = LN 2 ^ s+2 \ in contrast to L TF = L{NgL) l ^ s+1 ^ in the TF 
case. Condition (1.12) means in this region that Nr/L^ ~ N s /( s + 2 )r / L — > 0. Since 
-/W/Lll ~ (p/g)(a/r), this condition is automatically fulfilled if g/p is bounded away 
from zero and a/r — > 0. Conversely, if g/p is bounded, then (1.12) implies a/r — > 0. 
The energy functional is 

C>] = / (Vl(z) P (z) + p{zfe(g/p{z))) dz (2.19) 

JR 

with corresponding energy 

E LL (N,L,g) = inf {^[p] : p > , jf p(*)d* = ^} • (2-20) 
Introducing the density parameter 



7 := N/L LL = (N/L)N- 2 ^ S+ V 



(2.21) 
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we can write the scaling relation of the functional as 

E LL (N,L,g)=N 1 2 E LL (l,l,g/ 1 ) 
The minimizer, p LL (z), satisfies 

P^ L , g {l^z)= 1P \\ gh {z) . 



(2.22) 



(2.23) 



Theorem 2.4 (LL limit). Suppose r/L — > and a/r — ► 0, wi/t 5/7 > /raed as 
A" — > 00. T/ien 



(AT7 2 )- 1 [£^ M (JV,L,a,r)-JVe- L /r 2 ] -► £ LL (1, 1,0/7) 



(2.24) 
(2.25) 



and 

weakly in L 1 
2.5 The Girardeau- Tonks Region 

This region corresponds to impenetrable particles, i.e, the limiting case g/p — > 00 
and hence the formula (1.6) for the energy density. As in Region 4, the mean density 
is here p ~ 7 = (N/ L)N~ 2 ^ s+2 \ The energy functional is 



ff T W= / fa(*)p(*) + (tt 2 /3)pW 3 ) <** 



with energy 



E GT (N,L) =inf jff T [p] : p>0, / p{z)dz = N 



which can be written 



E GT (N,L) = AV-E Ui (l,l 



2 771GT/ 



The minimizer p G J(z) has the form 



(2.26) 

(2.27) 
(2.28) 

(2.29) 



pgrco=*-v T -n*)ii /2 , 

with // GT determined by the normalization. Note that its shape is different from that 
of (2.16), which makes it possible to distinguish experimentally the TF regime from 
the GT regime. The scaling relation for the minimizer is 

p GT L (L LL z)= 7 p G J(z) . (2.30) 

Theorem 2.5 (GT limit). Suppose r/L — > and a/r — > 0, with g/j — > 00 as 
AT -> 00. T/ien 



and 

weakly in L 1 (IR). 



(AT7 2 )- 1 [^ M (AT, L, a, r) - A^/r 2 ] - £ GT (1, 1) 
I-'pZ^z) - p G J(z) 



(2.31) 
(2.32) 
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2.6 Limiting cases of the general energy functional 

As already stated, the proof of Theorem 1.1 in Section 4 consists in comparing 
the ground state energy of the many-body Hamiltonian (1.1) with the ground state 
energies of the functionals defined in Subsections 2.1-2.5 in the various parameter 
domains. To link these special cases to the functional (1.7) it then remains to show 
that the ground state energy of (1.7) coincides with that of the functionals in Sub- 
sections 2.1-2.5 in the appropriate asymptotic limits. The proof of this follows the 
same pattern as the derivation of the 3D TF limit from 3D GP theory in [27, 28] 
and we shall here only give explicit proofs for the functionals in Theorems 2.3 and 
2.4 as examples. The limit theorems for the density are derived from the energy 
convergence by variation of the external potential in a standard way (c.f., e.g. [28]). 

Proposition 2.1. If N -> oo, NgL -► oo, but gL TF /N = gL(NgL) 1 /( s+1 ') /N -► 0, 
then 



(N/L 2 )(NgL)°K' + V E(N,L,g) - £ TF (1,1,1) 



and 



weakly in L 1 



(N/L TF ) 1 piv,L, 5 (£tf^) -► pJJ A (z 



(2.33) 
(2.34) 



Proof. With p the scaled density given by (2.11) the energy functional (1.7) can be 
written 



£[p] = (N/L 2 )(NgL) s /^ / ((NgL)-^^ \d^p\z)\ 2 
+ V(z)p(z) +p(z) 3 N(gL TF )- 1 e(gL TF /Np(z))yz . 
Now ie(l/t) < \ for all t [22], so 

~ (N/L 2 )(NgL) s /( s+lS> ~\ _1 £[p] 

< I {{NgL)-^ +2 y^\dy/W)\ 2 + V(z)p(z) + \~p(z) 2 ) dz 

Let 

3£ (z) = (2e)~ 1 eM-\A/e) ■ 
Define p = pj\ 1 * j £ . Then, since \dj £ \ = e~ l j £ and f p(z)dz = 1, 



(2.35) 



/ 



\d^fp\z)\ 2 dz < l/{4e 2 ) < oo , 

and (2.36) implies 

hmsup (N/L 2 ){NgLyl^ " E(N , L,g) < £?J\p% A * j e ] 



>n -l 



(2.36) 
(2.37) 

(2.38) 
(2.39) 
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in the limit considered. Moreover, f (pJJ : i * je) 2 < / (PiJ,i) 2 because f j £ = 1, and 
/ P$,i*Mz)\z\ a dz - I P ^z)\z\"dz (2.40) 



for e —> 0. (Note that Pi^i = [/« — M s ]+ is continuous and of compact support.) 
Hence 

limsup \(N/L 2 )(NgL) 8 K a+1> >]~ 1 E(N,L,g) < E TF (1, 1, 1) . (2.41) 
On the other hand, dropping the positive gradient term in (2.35) gives 

~ (N/L 2 )(NgL) s l( s+v> X~ 1 £[p] > f (V(z)p(z) + p{zf Me(l/M p{z))) dz , (2.42) 

with M = N/ (gL^p). Note that M — > oo in the limit considered here. The functional 
on the right side of (2.42) has as minimizer 

p iM \z) = [f M 1 (pW-V(z))] + (2.43) 

where f^ 1 is the inverse of the function /m(0 = d/dt[Mt 3 e(l/tM)} and ^ M ^> is 
chosen so that J p( M ^ = 1. Note also that i _1 /m(£) 1 as M -> oo, uniformly on 
[5, oo[ for every <5 > 0. From this it follows easily that converges uniformly to 
pJJi, given by (2.16), as M — > oo. With p = pN,L,g, the minimizer of £ , we thus 
obtain from (2.42) 



lim inf 

iV^oo 



(N/L 2 )(NgL) s ^ s+ ^ E(N,L,g)> E lt (1,1,1) . (2.44) 



-i 



TF/ 



To prove the corresponding result (2.34) for the density we pick a C°° function 
Y of compact support together with an e > and replace Vl{z) by 

V L (z) + ^(iV 5 L) s /( s+1 )y(L- 1 (iV 5 L)- 1 /( s+1 )z) = -^(NgLy/^ [V(z') + eY(z')] 

(2.45) 

with z' = z/L TF = L _1 (Ar 5 L) _1 /(s+i) z . While V(z') + £Y(z') is not strictly homo- 
geneous of order s, it is asymptotically homogeneous in the sense of Def. 1.1 in [28] 
and as in the proof of Lemma 2.3 in [28] this is sufficient for (2.33), now with the 
modified external potential (2.45). Since both (1.8) and the TF energy are concave 
in e, the derivative with respect to e can be exchanged with the limits N — > oo, 
(NgL) — > oo, giving (2.34) in the sense of distributions. Since the densities have 
norm 1, the convergence holds also weakly in L X (R). □ 

Proposition 2.2. If N -> oo with g/j fixed, where 7 = N/L LL = [N / L)N- 2 ^ s+2 \ 
then 

(N^E^^g) ^E LL (1, l, 5 / 7 ) (2.46) 

and 

7~W, 9 (W) - pi5 ifl/7 (*) (2.47) 

weakly in L x ( 
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Proof. With Lll = LN 2 ^ S+2 ^ we define the scaled density p by 

p(z) = (N/L LL )p(z/L LL ) . (2.48) 
The energy functional (1.7) can then be written 



£[p] = N 7 2 / N' 2 \d^/pXz)\ 2 + V(z)p(z) + p(zfe{g/( 1 p(z))) dz . (2.49) 
The lower bound 



(N^y'EiN, L, g) > E LL (1, 1, g/j) (2.50) 
follows simply by dropping the positive gradient term from the right side of (2.49). 

,g/i 



For the upper bound take j £ as in (2.37) and define p = j £ * p^\ , to obtain 



AT- 2 

(N^r'SiplK—^ + S^lp] (2.51) 



4e 2 

and hence 



iV-+oo 

for all e > 0. The convergence 



hmsup(iV 7 2 )- 1 J E(iV,L,< 7 ) < £\l gh [je *p\\ gh \ (2-52) 



lim /i L W7^ * * )9/7 ] = £i L w>i,i, 9/7 ] = ^ LL (1, 1,5/7) (2-53) 



follows by continuity of \z\ s and t e(t) and uniform convergence of j £ * PYig^ to 
LL 

The convergence of the densities follows as in the previous proposition by per- 
turbing the external potential, this time replacing Vl{z) by 

V L (z) + ej 2 Y (L -1 N~ 2 ^ s+2 ^ z) = 7 2 [V(z') + eY(z')] (2.54) 

with z' = z/L LL = L^N'^^z. □ 

2.7 One-dimensional GP as limit of three-dimensional GP 

We shall now demonstrate that the ground state energy in Regions 1-3 can be ob- 
tained as a limit of the three-dimensional Gross-Pitaevskii energy. The latter is 
defined by the energy functional 

4d>] = J (|V<I>(x)| 2 + {^(x- 1 ) + V L (z)] |<K(x)| 2 + 4H$(x)| 4 ) d 3 x . (2.55) 

We denote its ground state energy, i.e, the infimum over all $ with f |<1>| 2 = N, by 
E^(N,L,r,a). It satisfies the scaling relation 

E^(N, L, r, a) = (N/L 2 )E^(1, 1, r/L, Na/L) . (2.56) 

Because of (2.56) and (2.6) it is sufficient to compare E^r and E GP for AT = 1 and 
L = 1. 
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Theorem 2.6. Let g be given by (1.10). In the limit r — ► and a — ► 0, 

E 3 G D p (l,l,r,a)- e > 2 
# GP (M,<7) 

uniformly in g as long as r 2 E (1, l,g) — ► 0. 



(2.57) 



Proof. We denote the minimizer of the one-dimensional GP functional (2.4) with 
iV = 1, L = 1 and g fixed by <fi(z) 2 . Taking b r (x. ± )(/)(z) as trial function for the 3D 
functional (2.55) and using the definition (1.10) of g we obtain without further ado 
the upper bound 

E^(l,l,r,a)<e ± /r 2 + E GP (l,l,g) (2.58) 
for all r and a. For a lower bound we consider the one-particle Hamiltonian 

H rA = -A ± + V r ± (x ± ) - d\ + K(z) + 87ra6 r .(x ± ) 2 0(z) 2 . (2.59) 

Taking the 3D Gross-Pitaevskii minimizer <3?(x) for N = 1, L = 1, as trial function 
we get 

r (.2,2^2 



inf spec H r>a < E$g (l,l,r,a) - ina J <5> 4 + 8na J b 2 <p 2 <S> 2 
< E^(l,l,r,a) +4na J bffi 

= E GP (l,l,r,a) + 9 - f ^ . (2.60) 

On the other hand, inf spec£f rja can be bounded below by Temple's inequality [40], 
which says that for any Hamiltonian H with lowest eigenvalues Eq < E\ and expec- 
tation value (H) < E\ in some state, 

, , ((H-(H)) 2 ) , s 

Eo > (H) - U ^ ; . (2.61) 

We apply this to H = H r ^ a and the state defined by b r <p. Here 

{H) = t I +E GP (l,l,g) + 9 - J 4>\ (2.62) 

and since E\ > e^/r 2 with > e ± , (2.62) is smaller than E\ for r^E (1,1, g) 
small enough. Moreover, 

(H - (H))(b r <f>) = (8to6 2 - g)(t?b r (2.63) 

and thus 



((H - (H)) 2 ) = (/ ^ [(Svra) 2 (| b^ - , 2 ] < (| 6 ) (8™) 2 (| tf) ||6 r ||L 
< const. gUf^g j ^ < const. E GP (1, 1, g) 2 , (2.64) 
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where we used [28, Lemma 2.1] to bound ffH^H^o by the GP energy. Combining 
(2.60), (2.61) and (2.62) we thus get 

£ 3 G D P (1, 1, r, a) - e L /r 2 > E GP (1, 1, g) (l - const. r 2 E GP (l, 1, g)) (2.65) 

and the proof is complete. □ 

Remark. In combination with Theorem 2.2 this results demonstrates a fortiori that 
the three-dimensional GP limit theorem in [27] holds uniformly in r/L, provided 
NaL/r 2 stays bounded. A more direct proof of this fact, closer to the lines of 
[27], is certainly possible, but it requires redoing all estimates keeping track of the 
dependence on r/L. 



3 Finite n bounds 

Before we give the proof of our main Theorem 1.1 in Section 4, we will explain briefly 
the strategy, and give some auxiliary results in this section. In particular, we will 
derive upper and lower bounds on the ground state energy of (1.1) with the external 
potential Vl(z) replaced by a box with Dirichlet and Neumann boundary conditions, 
respectively. To obtain bounds on the full Hamiltonian (1.1), space will be divided 
in z-direction into small boxes of side length £, and the bounds of this section will 
be used in every box. The reason for this is twofold: first, this allows to consider 
an essentially homogeneous systems, without the additional difficulty of the external 
potential Vl(z), and secondly, by varying £ the particle number in each box can be 
controlled. This is necessary, since the bounds we obtain in every box will not be 
uniform in the particle number. 

Since the particle number in the boxes will be small (compared to N), we denote 
it by n. In this section, we study the Hamiltonian 

H = Y,(-^ + K ± (xf))+ £ ««(!*- x,-|) (3.1) 

j=i i<*<j<™ 

on L 2 ((R 2 x [0, i]) n ). Let E§ M (n,£,r,a) and E® M (n,£,r,a) denote its ground state 
energy with Dirichlet and Neumann boundary conditions, respectively. The following 
Theorem gives upper and lower bounds on the ground state energy of (3.1) in terms 
of the ground state energy of the ID Hamiltonian (1.3). Its proof will be given in 
Subsections 3.1 and 3.2. A crucial step in the proof of the lower bound will be the use 
of the 'Dyson Lemma' (see [9, 31]), which converts the 'hard' interaction potential 
v a into a 'soft' one. With this new 'soft' potential perturbation theory can be done, 
and rigorous bounds are obtained with the use of Temple's inequality (2.61). 

Theorem 3.1 (Finite n bounds). Let E^P (n,£,g) and E^p(n,£,g) denote the 
ground state energy of (1.3) on L 2 ([0,£] n ), with Dirichlet and Neumann boundary 
conditions, respectively, and let g be given by (1.10). Then there is a finite number 
C > such that 

'o\i/8l 



ET(n,e,r,a) - ^ > Bg> Ml ) (l - Cn (^f [l + VL (5) 



(3.2) 
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Moreover, 



E§ M (n,£,r,a) 



ne 



„2 — 



<E™(n,£,g) 



l + C 



at 

1 + ^2 



1/3^ 



(3.3) 



provided the term in square brackets is less than 1. 



Let us comment briefly on the error terms in (3.2) and (3.3). As already men- 
tioned, in the proof of Theorem 1.1 we will divide space in the z-direction into small 
boxes of side length I. The number of particles in each box will be roughly n ~ N£/L, 
where L = Nf p the extension of the system in z-direction. The n-dependence of the 
error term in (3.2) restricts us essentially to have a finite particle number n, i.e., that 
n ~ Nt/L ~ 1, or £ ~ L/N. But for (3.2) to be useful we need £ > r, i.e., r < L/N, 
or, in other words, r should be of the order of the mean particle spacing, or smaller. 
For r S> L/N, r is much bigger than the mean particle spacing, and we have to use 
a different strategy, similar to the one used in the 3D problem [31, 27]. This will be 
necessary for the lower bound in Regions 1-3. The result is stated in Theorem 3.2. 
For it's proof it will be necessary to use the box method also in x^-direction, similar 
to the 3D case considered in [27]. However, one cannot use directly the results from 
there, one has to be careful to retain uniformity in r/L. 

Likewise, (3.3) will not be useful for an upper bound in all the Regions 1-5. The 
reason is the last term in (3.3), where we want gl ~ gL/N < 1, which is only fulfilled 
in Regions 1-4. For Region 5, we use a different upper bound, given in the following 
Theorem. The proof of Theorem 3.2 will be given in Subsections 3.3 and 3.4. 



Theorem 3.2 (Additional energy bounds). 

Thru. 3.1, 



With the same notations as in 



E$ M (n,£,r,a)- 



ne 



n 2 g 
2£ 



1 - C 



n 



-1/14 



> 



+ 



nlaV'* 
r 2 J 



1 + 



Jnt 



na 



l/4\ 4 /39 



na 
+ T 



Moreover, denoting the range of v a by Rq, 

ne- 1 . vr 2 n 3 (1 + 1/n) (1 + l/2n) 
3~^ (i _ ( n _ l)R Q /tf 



E% M {n,l,r,a) 



< 



(3.4) 



(3.5) 



provided (n — l)Ro < £. 

Remark. By definition, Rq ~ a. Eq. (3.4) will be used with £ ~ (^L/N) 1 / 3 
and n ~ N£/L. In this case we have, in Regions 1-3, n£g ~ a(N/r 2 L) 1 / 3 <ti 1, 
na/£ ~ ap < 1, y/na/r ~ a(A^/r 2 Z) 1 / 3 < 1 and r/(^/n£) ~ 1. 

The following four subsections contain the proofs of Theorems 3.1 and 3.2. 
Throughout, C denotes a constant independent of the parameters, although the 
value of different C's may be different. 
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3.1 Upper bound for Theorem 3.1 

In this subsection we are going to prove (3.3). We use the variational principle. Let 
ip denote the ground state of (1.3) with Dirichlet boundary conditions, normalized 
by (i/)\if>) = 1, and let G and F be given by 

n 

G(xi,...,x n ) =ip(z 1 ,...,z n )Y{b r {x-f) (3.6) 

i=i 

and 

F(xi,...,x n ) = JJ/(|x i -x J -|) . (3.7) 

i<j 

Here / is a monotone increasing function, with < / < 1 and f(t) = 1 for t > R for 
some R> Rq. For t < Rwe shall choose f(t) = fo{i)/ fo(R), where /o is the solution 
to the zero-energy scattering equation for v a [31, 27]. Note that fo(R) = 1 — a/R for 
R > Rq, and /o(£) < min{l, a/t}. We use as a trial wave function for (3.1) the 
function 

*(xi, . . . ,x n ) = G(xi, . . . ,x n )F(xi, . . . ,x n ) . (3.8) 

Let p^ denote the two-particle density of ip, normalized by J p^ (z, z')dzdz' = 1. 
Since F is 1 whenever no pair of particles is closer together than a distance R, we 
can estimate the norm of by 

<*!*>> Jg 2 ^l-^0(i?-|x i -x j |)j 

J p^\z,z')b r (^)\(y ± fe(R-\^-y\)dzdz'd^ ± d 2 y 1 - 
J pJ ) (z,z')6,(x ± ) 2 6,(y ± ) 2 e(i?-|x ± -y ± |)^'^dV 
1 6 r (x^) 2 6 r (y^) 2 0(i?-|x^-y^|)d 2 x^V 

>1- ^^116111, (3-9) 

where we used Young's inequality [24] in the last step. 
Using 

(tf | - Aj\V) = - J F 2 GAjG + J G 2 \VjF\ 2 (3.10) 

and the Schrodinger equation H n ^ g ip = E^tp, we can write the expectation value of 
(3.1) as 

(mm = (< + ^)(^)-^i^^-^)i*) 



= i 
> i 
= i 



MS 



+ g 2 Ei v ^i 2 +E^- x ^i)i F i 2 • ^ 



i<j 
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Now, since < / < 1 and /' > by assumption, F 2 < /(|xj — x,,-|) 2 , and 

£|V,F| 2 <2^/'(|x l -x,|) 2 + 4 / / (|x fe -x i |)/'(|x fe -x J |). (3.12) 

j=l i<j k<i<j 

Consider the first term on the right side of (3.12), together with the last term in 
(3.11). These terms are bounded above by 

2E/ G2 (/ , (|x,-x J |) 2 + It; a (|x J -x J |)/(|x l -x,|) 2 ) 

i<j 

= n(n -I) J b r (^fb r (y^fpf{z, z>) (/'(|x - y|) 2 + \v a {\* - y|)/(|x - y|) 2 ) . 

(3.13) 

h(z) = J (/'(|x|) 2 + It, a (|x|)/(|x|) 2 ) d^ . (3.14) 

Note that h{z) = for \z\ > R, and J h{z)dz = 4ira(l — a/R)^ 1 . Using Young's 
inequality for the integration over the _L- variables, we get 

(3.13) < ni - n ~ l) \\b\\l [ pf(z, z')h{z - z')dzdz' . (3.15) 

Consider now the contribution from the last term in (3.12). We can write 

4 £ J G 2 /'(|x fe - XiD/'dxfc - x,|) = 2 -n{n - l)(n - 2) 

k<i<j 

■ J / / (|x 1 -x 2 |)/ / (|x 2 -x 3 |)6 r (x^) 2 6 r (x^) 2 6 r (x^) 2 pJ !) (z 1 ,z 2 ,Z3)rf 3 xid 3 x 2 d 3 X3 , 

(3.16) 



Let 



(3) 

where p^ denotes the three-particle density of ip, normalized by 1. Let 

k(z) = f /'(IxD^x 1 - , (3.17) 
which is supported in [-R, R]. For the integration over x^ we use 

J /'(|X! - x 2 |)6 r (x^) 2 dV < ^-Hz! - Z 2 ) < ^||fc||oo • (3-18) 
For the remaining integrations, we proceed as in (3.15) to obtain 

(3.16) < |n(n-l)(n-2)l^&M||Ar|| 00 / pf(z, z')k(z - z')dzdz' . (3.19) 
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Now, for any G if^R), 
\W(z)\ 2 -W)\ 2 \ 



dz' 



< 



f 

2\z-z>r([ \<P\ 2 ) 1/4 ([ 

\JR / \JR 



z>) 



dz' 

2\ 3 / 4 



(3.20) 



where we used H^H^ < ||(i0/(iz||2 ||</>||2- Applying this to p$(z,z'), considered as a 
function of z only, and using the fact that the support of h is contained in [-R, R], 
we therefore get 



f pf{z,z')h{z-z')dzdz'- f h{z)dz f pf(z,z)dz 



< 2R 1 / 2 [ h(z)dz [ 
Jr J 



dzdz' < 2R 1 ' 2 



[ h(z)d 
Jr 



z ( t/j 



' dz\ 



3/4 



(3.21) 



(2) 

where we used Schwarz' inequality, the normalization of p , and the symmetry of 
ip. The same argument is used for (3.19) with h replaced by k. Now J h(z)dz = 
47ra(l - a/R)~ l , and 



|oo < 



2na 



(1+lnGR/a)) , 



l-a/R 

/",.., 2iraR /_ a \ 



(3.22) 



Therefore 



(3.15) + (3.19)<^^^ 1 + K 



l-a/R 



j pf{z,z)dz + 2R 1 / 2 U 



dz\ 



3/4' 



where we denoted 



2vr (n _ 2) al + HRM (R 
3 1 j i? 1-a/i? 



loo • 



(3.23) 



(3.24) 



It remains to bound the second term in (3.11). We use again the fact that F is 
equal to 1 as long as the particles are not within a distance R. We obtain 



i<j 



> 



n(n-l) f ( 2 ) 



> 



2 

nln — 



| pf(z,z)dz J b r {^fb r (y^ 



l-6(R-\x A 



(3.25) 
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Putting together the bounds (3.9), (3.23) and (3.25), and using the fact that 
g\n(n — 1) J p^\z,z)dz < and (tp\ - d 2 /dz 2 \tp) < E^P /n, we obtain the upper 



bound 



-J^V) — D {n,e,9) 



l-^^^\\b\\i 



1 a/R + K ^R 2 imi4 



provided the term in square brackets is positive. We now use E^ /n > (ir/£) 2 , and 
choose 

(3.27) 



R? 



ar 



n 2 {l+gi) ' 

This gives (under the assumption (na/r) 2 (l + gi) < 1) 

E§ M (n,£,r,a) - ^ < <(M,<?) (l + C (^) 2/3 (1 +^) 1 / 3 ) (3.28) 
for some constant C > 0. 
3.2 Lower bound for Theorem 3.1 

We are now going to prove (3.2). We write a general wave function ^ as 



*(xi,...,x n ) =/(xi,...,x n )fj6 r (x^) , 



(3.29) 



k=i 



which can always be done, since b r is a strictly positive function. Partial integration 
gives 



_L n r 

<W> = ^r + £/ 
i=l J 



\V t f\ 2 + ±Y,Va(\^-Xj\)\ff 



k=l 



(3.30) 

Choose some R > Ro, fix i and Xj, j / i, and consider the Voronoi cell Qj around 
particle j, i.e., 



Qj = {x : |x — Xj\ < |x — x/j| for all A; 7^ j}. 
Denote by Bj the ball of radius R around Xj. We can estimate 



(3.31) 



/ fe r (x^) 2 (|V,/| 2 + l Wa (|xi - x,|)|/| 2 ) d 3 x, 
^minb^x^a / U(\^ - x,|)|/| 2 



_L\2 



> 



min xeBj . 6 r (x 
max x6B 6 r (x ± ) 2 



a/ ft^x^^f/dx.-x.DI/l 2 , (3.32) 
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where we used Lemma 1 of [31], and 

U(r) = { 3{R3 ~ ^ f ° r R ° ~ r ~ R (3.33) 
\ otherwise . 

For some 8 > define B$ C R 2 by 

B s = {x^ G R 2 : b^) 2 > 5} . (3.34) 

Estimating max xe ,g. 6 r (x J -) 2 < minxg^ 6 r (x J -) 2 + 2(i?/r 3 )||V6 2 || 0O , we obtain 

min xeB . b r (^) 2 ± / i? ||Vb2||^\ 

> XB 4 (x 7 - /r) 1 - 2- . (3.35) 



max xeB ^(x- 1 ) 2 



^3 



(For a proof that V6 2 is a bounded function, see the proof of Lemma 1 in the 
Appendix). Here xb s denotes the characteristic function of Bs- Denoting k(i) the 
nearest neighbor to particle i, we conclude that, for < e < 1, 



l[b r ( X i) 2 d 3 Xk 



k=l 



E/ |V,/| 2 + i ^^(Ix.-x.DI/p 

n ,. 

>E/ [e|Vi/| 2 + (l-e)|V,/| 2 Xm in fcN -z k \>R\Zi) 

n 

+a>U(\x i - x fc(i) |)xB 4 (x£ (i) /r)|/| 2 ] J] b r (^) 2 d^ k , (3.36) 



i=l 



k=l 



where a' = o(l - e)(l - 2i?||V6 2 || 00 /r<5). 
Define . . . , z n ) > by 

/n 
|/(x 1 ,...,x n )| 2 n^(x^)Vx^ • (3-37) 
k=l 

Neglecting the kinetic energy in _L-direction in the second term in (3.36) and using 
the Schwarz inequality to bound the longitudinal kinetic energy of / by the one of 
F, we get the estimate 

(W)-^> 

n „ n 

E / [e|^| 2 + (l- e )|^F| 2 XminfcN -z k \>R(zi) \ dz k 

i=l J k=l 

n „ n 

+ E / [elV^/l^a'C/d^-x^Dxs^x^/r)!/! 2 ] Hbr^ifd^ , 
i=i ^ fe=i 

(3.38) 
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where dj = d/dzj, and V -1 denotes the gradient in _L-direction. We now investigate 
the last term in (3.38). Consider, for fixed zi, . . . , z n , the expression 

n „ n 

£ / [ £ I V ^! 2 + a ' U ^ ~ x«i)\)XB,{x£ (i) /r)\f\ 2 ] J] bM) 2 d 2 *i • (3.39) 
i=l J k=l 

To estimate this term from below, we use Temple's inequality, as in [31]. Let e 1 - 
denote the gap above zero in the spectrum of — A 1 - + V 1 - — , i.e., the lowest non-zero 
eigenvalue. By scaling, e L /r 2 is the gap in the spectrum of —A 1 - + V^h — e^/r 2 . Note 
that under the transformation § ^b~ x § this latter operator is unitarily equivalent 
to V^* • V- 1 as an operator on L 2 (R 2 , 6 r (x ± ) 2 d 2 x ± ), as considered in (3.39). Hence 
also this operator has e lL /r 2 as its energy gap. Denote 

(TL \ ^ TL 

E^lx-x^DxB^x^/r) [[Mx^^- (3-40) 
i=i / fc=i 

Temple's inequality (2.61) implies (under the assumption that the denominator in 
the last term is positive) 

(3.39) > |F|V <t /> (l - ^ g[ - a L_ ) . (34!) 

Now, using (3.33) and Schwarz' inequality, (U 2 ) < 3n(R 3 - R^iU), and 

(U) < n(n-l)Ju(\x-y\)W(x ± ) 2 b r (y L ) 2 cPx J -<1 l y L 

< n(n - 1)M J t/dx,)^ < n{n _ 1)^ . (3.42) 

Using this and a' < a in the error term, we obtain 

(3.41) > \F\ 2 a"(U) , (3.43) 

where 



a" = a' 1 



3n ar 2 1 



l-^-"37r||6|li 



-IN 



(3.44) 



ee- 1 i? 3 1 - (iV-ft) 3 L ee X R " 114 1 - (R /R) 3 _ 
with the understanding that the term in square brackets is positive. Now let 

d(z - z>) = f 6 r (x ± ) 2 6 r (y ± ) 2 t/(|x - y|)^(y ± /r)d 2 x ± d 2 y J - . (3.45) 

Note that d(z) = if \z\ > R. We estimate (U) from below by 

„ n 
i^j fe, ky^i,j 1=1 

> ^Ju^-^DxB^f/r) (l-E ^(^-|x fe -x J |))n6 r (x i -) 2 ( iV 

> Y J d{z i -z j )U-{n-2)^\\b\\l\ . 



(3.46) 
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Let 



7r/? 2 

o M -n"[l-( n -2)^-\\b\L, 



(3.47) 



and denote g' = 2a!" J R d(z)dz. Note that since ^(x- 1 ) 2 - 6(y ± ) 2 | < -R||V6 2 ||oo for 
|x ± -y ± |<i?, 



f d(z)dz > %(\ b{^fd 2 ^ - i?||V6 2 ||ooA 



47T 



> (\\ b \\i-5-R\\Vb 2 \Ur) 



We write 



E 



e\d i F\ 2 + a'"Y J d{z i -z j )\Ft 

/-l It 
-^ T \d t F\ 2 + a"'d(z i -z j )\F\ 2 J[dz k 

i^j L k=l 

Now consider, for fixed Zj, j / i, the expression 

/ ^\d i F\ 2 + a. m d(z i -z j )\F\< 



fc=i 

n 



We claim that 



(3.50) >y max |F| 2 X[iV _ R] (^) 1 
\zi-Zj\<R \ 



Assume that (3.5f) is false. Estimating, for any i? 1 ([0, £])-function 0, 



\Hz)\ 2 - max |0(z)| 2 > 

\z-zo\<R 



d\4 



> -2R 1/2 max \6(z) 
\z-z \<R 



(3.48) 



(3.49) 



(3.50) 



(3.51) 



1/2 



(3.52) 



and applying this estimate to F (considered only as a function of Zi), we obtain, 
using ef\diF\ 2 < ±(n - l)g' max\ z ._ z .\< R \F\ 2 by assumption, 

a'" J d(z t -z 3 )\F\ 2 dz t 

> a'" / - z 3 )dz v max ^ |F| 2 (l - - itf *) ^ 



> 5SX[H,<-fl]( z J'), max |F| 2 1 



n 



\ i/2 N 



(3.53) 



contradicting our assumption. This proves (3.51). 
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Putting everything together, we thus obtain 

j_ n n 

{mm - ^ E / [(! - e )l 5 ^| 2 Xmin fc |*-* fc |>*(*i)] II d ** 

i=l ^ fc=l 

+ /, max iFpx^i?]^) ]J ' ( 3 - 54 ) 

where 

5" = 5' ^-^(n-iyi?) 1 ^ . (3.55) 

Assume that (n + < £. Given an F with J |F| 2 dzi • • • dz n = 1, define, for 
< zi < z 2 < ■ ■ ■ < z n < £ ~ (n + 1)R, 

i/j(z!, ...,z n ) = F(z! + R,z 2 + 2R, z 3 + 3R,...,z n + nR) , (3.56) 

and extend the function to all of [0,1— (n + l)R] n by symmetry. A simple calculation 
shows that, for 

n 

H' = (l-e)J2 -9? + 9" £ S( Zi - Z j) (3.57) 

i=l i<j 

on L 2 ([0, £ — (n + l)R] n ), 

(3.54) > WH'W) > (l-e)E™(n,£-(n + l)R,g")(m 

> {l-e)E™{nA9"){m • (3-58) 

Here (ip\H'\tp) is interpreted in the quadratic from sense, since if) does not necessarily 
fulfill Neumann boundary conditions. Since these give the lowest energy for the 
quadratic form, (3.58) is valid anyway. 

It remains to estimate (if>\ip) for the F that is related to the true ground state ^ 
by (3.37). We have 



„ n 

<</#} = J \F\ 2 Hxm-R](zk)lle(\ z i- z j\- R ) 



I 



k=l i<j 
n 



> 1 - / \F\ 2 



Yj i 1 - X[R,e-R](.Zk)) +J2 6 ( R - \ Z i - z i\)X[R,t-R\{Zj) 
k=l i<j 

(3.59) 



To bound the second term on the right side of (3.59), we use 

Y I \ F \ 2d (R~ \* ~ Zj\)X[R,e-R](zj) < RY, [ , ma ^J F \ 2 X[R,t-R](z 3 ) 

< y(Er(N,£,r,a)- ^f) , (3.60) 
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where the last inequality follows from (3.54). For the other term, we use the simple 
fact that, for any function <j> 6 i? 1 ([0, £]), and for < R < £, 

f R ^ z)dz _ * [ ^ Z ) dz = ( L <j>>(z)f R/ (z)dz , (3.61) 
Jo 1 Jo Jo 

with fn/{z) = zR/t — min{z,i?}. Note that \ fn y e{z)\ < R- Applying this to F 2 and 
using Schwarz' inequality we get the estimate 

r n R /I \ 1/2 

J \F\*Y,{l-Xvu-q{*k)) < 2n^ + 4nR(-(F\Z t -d?\F)) 

R / 1 „_L\V 2 

< 2n^ + AnR[-E® M (n i e,r,a)-^ v ) . 
I \n r A ) 

Denoting A = E® M (n, £, r, a) — ne^/r 2 , we thus have 

A > E™(n,£,g") (l - ™A - 2n| - 4 v ^iL4 1 / 2 ) , 



(3.62) 



(3.63) 



which implies 



A > 



l + 2RE^{n,£,g")/g' 
We now use the simple upper bounds 



^l-2nj-2R^nE^(n,£,g")^ . (3.64) 



7r 2 n 3 



E^{nA,g")<^^g" and E^{n^)<~ 



(3.65) 



which follow from a constant and a free-fermion trial wave function, respectively. 
Moreover, by concavity of E^P in g, E^p{n,£,g") > E^{n,£,g)g" /g for g" < g, and 
therefore 



E^ l (n,£,r,a) 2 



P QM 

We now choose 



>E^(n,£,g) 9 -[ 1 



R 



^ 27T 2 

nin + 1) H i=n 

V3 



a\V8 



(3.66) 



(3.67) 



and obtain 



<V ^, a) - ^ > <(n, <, 5 ) (l - Cn (£) V§ [l + f $ ^ ) (3.68) 



for some constant C > 0. 



LSY August 14, 2003 



27 



3.3 Upper bound for the Girardeau- Tonks regime 

We use as a trial function 

n 

tf(xi,...,x n ) = il>(z 1 ,...,z n )Y[b r (xk) , (3.69) 

k=l 

where ip is the ground state of a one-dimensional Bose gas of particles with hard cores 
of radius i?o- It is well known that it's energy is the same as that of n non-interacting 
fermions on a line of reduced length t — (n — 1)Rq- An explicit calculation yields 
(3.5). 



3.4 Lower bound for Theorem 3.2 

We start by defining an auxiliary 2D GP energy functional by 

/" + V ± \<P\ 2 +p\ ( j ) fyd 2 ^ (3.70) 

Jw? v J 

for some parameter p > 0. The following fact will be needed below. 

Lemma 1. For any p > 0, there exists a unique minimizer (up to a constant phase 
factor) of (3.70) under the normalization condition f \4>\ 2 = 1, denoted by (f> p , that 
can be chosen strictly positive. Moreover, both 4> p and V^~(f) p are bounded uniformly 
in p and x^ for p in a finite interval [0, P] . 

The proof of this Lemma, as well as the proof of Lemmas 2 and 3 below, will be 
given in the Appendix. 

The energy corresponding to the minimizer (f> p will be denoted by E aux (p). Define 
also (j> P:r ('x. ± ) = r _1 p (x ± /r). Writing the wave function as 

n 

*(xi, . . . ,x n ) = /(xi,.. ,x n ) Y{ 4>pA*k) > ( 3 - 71 ) 

k=l 

using partial integration and the variational equation for (j) p , we obtain 



i=i 



|Vi/| 2 + \ "-(I* ~ ^'Dl/I 2 " 2p^ p>r (Xi L ) 2 |/| s 



]j0 p , r (x^) 2 d 3 x fc . 



fc=l 



(3.72) 



We now divide space in x^-direction into boxes of side length s, labeled by a. Let 
0a,max and <^ a ,min denote the maximal and minimal value of 4> Ptr inside box a, re- 
spectively. We obtain 

E$ M (n,e,r,a) - ^E^{p) - ™f / $ > inf T, [ E -M - 2n a p< max ] , (3.73) 
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with 

" L [l^/| 2 + lE J , J y^g(|x,-x J |)|/| 2 ] niiV(x^) 2 ^ 

(3.74) 

In (3.73) the infimum is over all possible distributions of the n particles into the 
boxes a. As a first step we will replace </> 2 max by 4> 2 a min in the last term in (3.73). 
The error in doing so is bounded above by 

2np sup (<f> 2 - ^ ) < 2v / 2ng||V ± ^|| 00 . (3.75) 

a ' 

To bound E a from below, we divide a into even smaller boxes, denoted by /3, 
with side length t < s, where s/t £ N. Let q denote the number of boxes with 
exactly i particles. Then 

E a (n) > inf V a inf Ep(i) . (3.76) 

Note that the infimum is under the constraints ^ q = (s/t) 2 and Yl c ^ = n - For a 
lower bound on Ep, we use 

W ' fcj § PFIA ' 

(3.77) 

Fix some 5 > 0, and assume that (f) 2 a min > <5/r 2 . Then 



/3,max / 



(3.78) 



(compare with (3.35)). The rest of (3.77) can be bounded below by the same method 
as in the previous subsection. The only difference lies in the fact that b r is replaced 
by the constant function, and the integrations are only over the box (5. The result is 
(compare with (3.2)) 

Ep(n) > E^ D (n,£,8-ira/t 2 ) 

•(l-On( ? ) [l + -( 7 ) jj^l-^LJ-J. (3.79) 

To proceed we need an explicit lower bound on that will be proved in the 
Appendix. 

Lemma 2. There is a finite number C > such that 

Eh D (n,l,9) > l^p^O (l - Cn(i g f' 2 ) . (3.80) 
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Applying this lemma to (3.79), we obtain 

,_, , , n(n -l)47ra/ ^ /axVsT nt/a\ i/8]\ 



. M - y^ni lZ^i j (l _ Cnffa/t 2 ) 1 ' 2 ) . (3.81) 

Note that the right side is independent of (3 C a. We insert this bound in (3.76), and 
use the following Lemma. It is a simple generalization of a result of [31]. Although 
we need at this point only the version proved in [31], we state the Lemma in this 
general form for later use. 

Lemma 3 ([31]). For n G NU {0}, let E(n) be a sequence of non-negative real 
numbers that is superadditive, i.e., E(n\ +71-2) > E(n\) + E(ri2), and bounded below 
by 

E(n) > L(n)K(n) , (3.82) 
with K : W + — ► W + monotone decreasing, L : M + — ► R + convex, L(0) = 0, and 

L\x) < ^ (3.83) 

for some A > 1 and all x > 0. (Here L' denotes the right derivative of L.) Let c n be 
a sequence of non-negative real numbers, with 

c n < M and ^ c n n = N . (3.84) 

n>0 n>0 

Then 

Yc n E{n) > ML(N/M)K((XN/M}) , (3.85) 

n>0 

where (x) denotes the smallest integer > x. 

The proof is given in the Appendix. Note that inlp Ep(n) is certainly a super- 
additive function, as the infimum over superadditive functions. Therefore we can 
apply Lemma 3 with L(x) = x[x — 1] + and A = 4, to (3.76), together with (3.81), 
and obtain 

4iran 2 / Is 2 ' 



where 



E a (n a ) > 11 - ) K(n a ) , (3.8(5) 



*(») = (l-C(4ntV S 2 >0 1/8 (l + <4ntV S 2 >^© 1/8 )) 

• (l - V2(Ant 2 /s 2 ) i l|V± ^ l|oo N j (1 _ C (4nt 2 /s 2 ) {ia/t 2 ) 1 ' 2 ) . (3.87) 
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We now insert this bound in (3.73), and use n a < n in the error terms. Note that 1Z 
is monotone decreasing in n, hence lZ(n a ) > lZ{n). Minimizing over n a gives 



E a (n a ) - 2n a p4> 2 in > 



^W. miri 2vrar 2 \ 2 1 



Ana 



1 + 



£t 2 p5 ) K(n) 



(3.88) 



This holds for boxes a where 4> 2 a min > 5/r 2 . In boxes where this is not the case, we 
simply use positivity of E a (which follows from positivity of v a ) in the form 

E a (n a ) - 2n a p(j) 2 a min > -2n a ^ . 

Putting everything together, using J2 a mm — r ~ 2 I ftp an< ^ choosing 

Airan 



P 



(3.89) 



(3.90) 



we obtain 



E® M (n,£,r,a) > ^E^iAiran/t) 



n 



Ai:an 2 I /-,,„ i , o ,. s 
~p- K8||V^|U- 



p i, 0O - + 2(5 + 



1 + 



2 \ 2 , 

r z \ 1 



2t 2 5nJ K(n) 



. (3.91) 



We are still free to choose the parameters t, s and 5. 

It remains to derive a lower bound on E anx . This will be done similarly to the 
lower bound on the 3D GP energy given in Subsect. 2.7. Consider the auxiliary 
Schrodinger operator 



ff aux = -A ± + V ± (x ± ) + 2pb{^f . 



(3.92) 



Using 4>„ as a trial function, we have 



infspec# aux < E anx (p) -pj $ + 2pJ b 2 <j) 2 p 
< E anx (p) +p Jb 4 . 



(3.93) 



On the other hand, using Temple's inequality (2.61), 



infspec# aux > e ± + 2p 



Ap 2 J b 6 



> + 2p / b" 1 



2pfb* 
2p\\b\\ 



2 

oo 



e±-2p\\b\\i c 



(3.94) 



Eqs. (3.93) and (3.94) together give 
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We now choose 



er, t = — — , 5 = e , (3.96) 



and 

l/4\ 4 /39 



1+ ' 



'no 



(3.97) 



and obtain as the final result 

E^(n,tr, a) >^ + ^(l-C[ e + ^}) . (3.98) 

Note that we did not pay any attention to the restriction s/t G N in choosing s and 
t. However, since, with our choice, s/t = e 2 y/n > n 1 / 2-1 / 7 , this can easily be made 
an integer by replacing e by some e with e < e < 2e. This affects only the constant 
C in (3.98). 

3.5 Boundary conditions 

As a last step in this section, before giving the proof of our main Theorem 1.1, 
we investigate the dependence of the ground state energy of (1.3) on the boundary 
conditions. In the upper bound above we used Dirichlet boundary conditions for the 
energy E 1T) , and Neumann boundary conditions for the lower bound. To relate these 
energies to the energy with periodic boundary conditions and to prove independence 
of boundary conditions in the thermodynamic limit, we need the following Lemma. 

Lemma 4. Denote Ep D (n,£,g) the ground state energy of (1.3) with periodic bound- 
ary conditions, i.e., on the torus [0,£] n . Then there is a finite number C > such 
that 

Eh D (n, £, g) < El D (n, I, g) < (n, i, g) , (3.99) 

and 

n 7/3 

Eh D (nJ,g) <E™(n,e,g) + C— • (3.100) 

Proof. The first inequality (3.99) is standard, noting that the interaction considered 
has zero range. For (3.100) we use a result of [37, Lemma 2.1.13 and Prop. 2.2.10], 
which implies, for < b < £/2, 

In 

E™(n,£ + 2b,g)<E™(n,£,g) + -p . (3.101) 

Using E^(n,£+ 2b, g) = E^(n,£,g(l + 2b/£))(l + 2b/£)- 2 > E^(n,£,g){l + 2b/£)- 2 
and E^(n,£,g) < ir 2 n 3 /(3£ 2 ), we obtain 

0,71 T)^ A-TT^ 1 h 

E™(n,£,g) < E™(n,£,g) + p- (1 + 2b/£) 2 + ^ — - (1 + b/£) . (3.102) 
Now b = const. £/n 2 ^ 3 gives the desired result. □ 
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4 Proof of Theorem 1.1 

With the results of the previous section in hand, we can now give the proof of our main 
Theorem 1.1. The proof will be divided into four subsections, two for the upper and 
lower bound, respectively. In each subsection, we compare the ground state energy 
of HN^L,r,a with the ground state energy of one of the functionals in Subsections 2.1- 
2.5, which, as explained there and in Subsection 2.6, is asymptotically equal to 
E(N,L,g) in the respective parameter region. Combining all the bounds obtained, 
this will prove Theorem 1.1, together with the claimed uniformity in the parameters. 
The corresponding convergence of the ground state particle density, as stated in 
Theorems 2.1-2.5, follows in a standard way by variation with respect to the external 
potential Vl{z) (compare with Props. 2.1 and 2.2 in Subsect. 2.6). Since the proof 
of the energy convergence is already quite lengthy by itself, the simple modifications 
necessary for a proof of the density convergence will be omitted. 

Let again L = N/ p denote the extension of the system in z-direction. As already 
mentioned in the beginning of Section 3, it will be necessary, for the lower bound to 
EQ M , to consider the case of small and large Nr/L separately. We will divide space in 
z-directions into small boxes, and use the bounds of Section 3 in each box. To control 
the number of particles in each box, Lemma 3 will be essential. For small Nr/L, 
where r is smaller than the mean particle distance, we will use the lower bound given 
in Thm. 3.1. For larger values of Nr/L, where r is actually bigger than the mean 
particle distance, the lower bound of Thm. 3.2 will be used instead. Note that this 
distinction is only relevant in Regions 1-3, since in Regions 4 and 5 Nr/L = rp <C 1 
by condition (1.12). Hence r is always smaller than the mean particle distance in 
Regions 4 and 5. 

4.1 Upper bound for Regions 1-3 

For an upper bound that gives the right asymptotics as long as g/p <C 1, we can 
essentially use the same technique as in [27] . The results of Section 3 are not needed 
in this case. As a trial function, we use 

N 

*(xi,...,xjv) = F(x 1 ,...,x JV )[]6 r (x^)0 GP (z A; ) , (4.1) 

fc=i 

where ^ GP = {pffig) 1 ^ 2 (cf. Subsect. 2.2), and F is the 'Dyson wave function', 
described in [9, 27]. The result is 

^ M (iV, r, a) < ^ + i? GP (iV, Z., 5 ) (l + Ca||6 r ||^ 3 ||0 GP ||^ 3 ) , (4.2) 

as long as a||6 r ||^ 3 ||</> GP ||TO 3 < 1. Note that, by the same proof as in Lemma 2.1 of 
[28], gU GF \\lo < 2E GP /N, and therefore 

E^ M (N, L,r,a)<^- + E GP (N, L, g) (l + const. a^N'^E^ (N, L, g) 1 / 3 ) . 

(4.3) 
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Now, in Regions 1-3, E GP (N,L,g) ~ E(N,L,g), and 



a 2 E GP (N, L, g)/N ~ a 2 (L~ 2 + S p) ~ (a/L) 2 + ^ (gpr 2 f < 1 . (4.4) 

p 



4.2 Upper bound for Regions 4 and 5 

If g/p is not small, the method of the previous subsection does not work, and we have 
to proceed differently. As in Subsect. 3.5, let E^(n,£,g) denote the ground state 
energy of (1.3) on an interval of length £, with periodic boundary conditions, and 
write E^ D (n,£,g) = n 3 e n (g£/n)/£ 2 . In [22] it is shown that, for every fixed t > 0, 
lim n _ >00 e„(t) = e(t). Since the functions are monotone increasing, concave and 
bounded in t, the convergence is actually uniform in t. By Lemma 4 of Subsection 
3.5, the same is true with instead of Ep D . Hence we get the estimate 

3 

^ D (n,£,g) < ^ (e{g£/n) + 8{n)) (4.5) 

for some bounded function 5 satisfying ]im n ^ 00 S(n) = 0. Without loss of generality 
we may assume that 5 is monotone decreasing. 

Let p be the minimizer of the LL functional (2.19) under the normalization con- 
dition J p = N. Note that p has compact support, with radius R = L(pL 2 ) l / s , 
where p = dE LL (N, L, g)/dN. (This R is different from that in Eq. (3.27).) By 
monotonicity and concavity of E LL in g, and by the scaling relation (2.22), 



S LL (iY,L,^ < ^ <3 i? LL (iV,L,( ? ) 



N 



N 



(4.6) 



Divide M in z-direction into intervals of length £, labeled by a, with Rq < £ < R/2. 
Let n a £ N be a collection of integers such that J2 a n a = Let V a = sup z£a V~l(z), 
and denote E® M (n,£,r,a) the ground state energy of (3.1) in a box of side length 
£ and with Dirichlet boundary conditions. By confining the particles into different 
boxes of length £ — Rq, a distance Rq apart, we get the estimate 



E QM (N, L, r,a)<J2 [E$ M (n a ,£ - Rq, r, a) + V a n c 



(4.7) 



Using (3.3) and (4.5) as well as monotonicity of e, we obtain 

£ 2 



E°- M (N, L,r,a)- — <Y^ -f (e{g£/n a ) + 5(n a ))n{n a ) + V a n a 



(4.8) 



with 



K(n) = 



(i - Ro/£) 2 



l + C 



/na\ 2 _ „. 

It) {l+9l) 



l/3^ 



(4.9) 



provided we choose n a and £ such that the term in square brackets is less than 1. 
Note the additional factor (1 — Rq/£)~ 2 , which is due to the fact that the size of the 
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box is only i — Rq. Now let n a = J p(z)dz, and n a = (n a ), where (x) denotes the 
smallest integer > x. With this choice J2 a n & — ^, but by monotonicity of (4.8) in 
iV we can plug in these values of n a for an upper bound. 

Since n a < n a + 1, and e and 5 are monotone increasing and decreasing, respec- 
tively, we obtain 



(4-8) < J2 



n « (e(gl/n a ) + «5(n Q )) + ft(*||p||oo + 1) + V a (n a + 1) 



(4.10) 

Here we estimated n a by ^||p||oo i n 72- • Denote V a = min z6Q , Vj, (z) . We estimate 
V a < V a + const. L- 2 (IIL){RIL) S ~ X in boxes where n a > 0. Using {R/L) s < 
3L 2 E LL (N,L,g)/N (see (4.6)), we therefore see that the error in replacing V a by 
V a is, in total, bounded above by const. E Lh (N, L, g)(£/R). 

Fix some < e < 1. We first bound the contribution to (4.10) from boxes 
where n a < 1/e. Now both e and S are bounded, the number of boxes with nonzero 
n a is bounded by (R/l) + 2, and V a < L~ 2 (R/L) S in these boxes. Therefore this 
contribution is bounded above by 

const. ^ (±R.{!\\pU + 1) + ^ LL( ^ L ' g) ) • (4-11) 

For the remaining boxes, we use n a > 1/e and n a < £\\p\\oo to obtain 

(4.10) < (4.11) + const. E LL (N,L,g) 4 

R 

- 3 

+ ( 1 + £ )-E ^( e (^/n a ) + 5(lA))(l + e) 2 ^llp||oo + l) + U a n a . 

a L J 

(4.12) 

Since x i— > x 3 e(l/x) is convex, we can use Jensen's inequality to bound the sum from 
above by 



/ 



P (zf (e(g/p{z)) + 5(1/ e)) (1 + ef ft(*||p|U + 1) + V(z)p(z)] dz . (4.13) 



Now, by the scaling (2.23), p(z) = jp g/l (jz/N), where 7 = (iV/L)iV" 2 /( s + 2 ) and 
pgl^ is a function that depends only on 51/7. Therefore ||p||oo < T 1 1 P^/t 1 1 00 5 and 

HpII! < ^7 2 IIp 9/7 IIL- 

We choose, for some < e < e, £ = (£7) 1 , and use 

^(1,1,3/7)^ <^<k 2/s ■ (4.14) 
Putting everything together, we get, for C denoting some universal constant, 

E^ M (N,L,r,a)- K ' A 



< (E^(N, L,g) + iV 7 2 (\\p gh \\lS(l/e) + C(676) 3 ) ) ■ K' , (4.15) 
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where 



K' =[\ + tt- —rr- + C- (1 + e) 




r 2 ^ 

The choice of e and e is determined by a/r, g/*y and N. The bound is uniform in 
g/j for bounded 5/7 and 7/5. 

If<7/7^ooasiV^oo (Region 5), we can use the same method to obtain an 
upper bound, replacing the bound (3.3) by (3.5) in (4.8). This gives a bound uniform 
in 5/7 (for 7/5 bounded). Combined with the result (4.15), this shows that in the 
limit TV — > 00 and r/L —> 



E^ M (N,L,r,a) - Ne L / 



2 



uniformly in the parameters, as long as a/r ^ and 7/5 stays bounded. 
4.3 Lower bound for Regions 3—5 

We now derive a lower bound on EQ M that will give the right asymptotics in Regions 
3-5. As in the upper bound, given in Subsection 4.2, we will use the box method, 
this time with Neumann boundary conditions. In each box, the results of Section 3 
will be used. In analogy to (4.5) we infer from [22] and Lemma 4 that 

3 

Eh D (n,£,g) > ^{e{gl/n) - S(n)) (4.18) 

for some bounded, monotone decreasing function 6 satisfying lim n >OQ S(n) = 0. This 

will be used in the bound for Regions 4 and 5. If gijn is small, however, we use 
e(g£/n) < \gtjn and (3.80) to obtain 

Eh D (nJ,g) > n2{n ~ l) e(gl/n)(l - Cn(gl) 1 / 2 ) . (4.19) 

We divide R in z-direction into intervals of side length M, labeled by a. Denote 
E® M (n, M, r, a) the ground state energy of (3.1) in a box of side length M and with 
Neumann boundary conditions. Let again V a = inf 2eQ Vl(z), and V a = sup zgQ Vl(z). 
By confining the particles into different boxes and neglecting the interaction between 
different boxes, we get the estimate 

E^ M (N, L, r, a) > inf £ [E$ M (N a , M, r, a) + V a N a ] , (4.20) 

a 

where the infimum is over all distributions of the N particles among the boxes a. As 
in the upper bound, we can estimate the difference of the maximum and minimum 
of Vl for boxes alpha a inside some interval [—R,R] by const. L~ 2 (M/L)(R/L) S ~ 1 . 
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For boxes outside [—R,R] we use V a > V a (l — sM/R). Choosing R the radius of 
the LL minimizer, and M = eR for some < e < 1, we see that, analogously to the 
upper bound, the error in replacing V a by V a (l — se) is, in total, bounded above by 
const.e£ LL (iV,L,5r). 

We now have to estimate E® M (N a , M, r, a) from below. We cannot directly use 
(3.2), because this bound is not uniform in the particle number. Instead we proceed 
similarly to [31]. We divide the interval M again into smaller intervals of length 
£ = rjM, where 1/rj £ N. Neglecting the interaction between different boxes, we 
obtain 

N a 

E$ M (N a ,M,r,a) > inf Y^CnE^inJ^^) , (4.21) 

where c n denotes the number of boxes containing exactly n particles, and the infimum 
is over all c n £ N under the condition ^ n c n n = N a and 2~^ n c n = M/£ = n^ 1 . 

Fix some < x < 1, and consider the case n > 1/x- Then, using (3.2) and (4.18), 



\n, £, r, a) - ^- > ^e(g£/n)(l - i/(n)) , (4.22) 
with 

nr /a\ 1 / 8 



, \ oil x) ^ /a\ i /* nr /a\ 
u(n) = ; [,\ + Cn - 1 + — ( - ) 
V ' e (at n) \rJ i \rJ 



(4.23) 



e(g£/n) \r> 

Note that i/(n) is monotone increasing in n. We now use Lemma 3 from Subsection 
3.4, with L(x) = x 3 e(g£/x). Note that for this L (3.83) holds with A = 6, since e is 
a monotone increasing and concave function, with e(0) = 0. Let N' = Yln>i/x CnU - 
The contribution from n > 1/x to the sum in (4.21^ will be bounded below using 
Lemma 3 and (4.22). For n < 1/x, we simply use (n,£,r,a) > ne L /r 2 . We thus 
obtain 



/3 



E$ M (N a ,M,r,a) - > —e{gM/N')(l - v((6N>n))) . (4.24) 

Using N a > N' > N a — 1/ (rjx) , this gives 

E$ M (N a ,M,r,a)-^ > ^e(gM/N a ) (l - ^J-)* (l-u((6N a n))) . (4.25) 

Now if eN a > 2 for some < e < x, we can choose |e < S < e such that 5N a G N, 
and take n = {SNa)' 1 . We also use {6/5) < (12/e) < 13/e (since e < 1). Note that, 
using (4.14) and 5 < e, 

^gt > ^^ LL (i, W7) 1/s = <t(gft) ■ (4.26) 

Therefore 

E$ M (N a ,M,r,a) - > ^(gM/N^ll , (4.27) 
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with 



(4.28) 



n = ( i - iy d - 5{l/x) - c - n 1/8 V + - f-V /8 ' 

V X/ V e(e^(flf/7)) e Vr/ [_ ^{g/l)r\rJ 

Here we used (4.14) and iV Q < iV in the last error term. 

The bound (4.27) holds for eN a > 2. If N a < 2/e, however, we use 

^e( 5 M/iV a ) < iV a ^ . (4.29) 
Using these bounds in (4.20), we obtain 

E^(N, L,r,a)-^ + CeE^(N, L, g) + 



> 



inf y 



N 3 

^ e (gM/N a ) + V a N a 



K(l - se) . (4.30) 



Note that, by (4.14) and (2.22), 



F = ^ s£ >' l »»WS)W' (4 - 31) 

Define p(z) = Y^ a N a x a (z), where Xa is the characteristic function of the interval a. 
The sum in (4.30) is bounded below by £ LL [p] > E LL (N, L, g), and therefore 

E^ (N , L , r , a) _ > E^(N,L, 9) K (l - C C - ^ E ,^ Xgh)Ml . ) ■ 

(4.32) 

The choice of e, e and \ is determined by 5/7 and a/r. They can be chosen such 
that (4.32) gives the correct lower bound in the limit considered, uniformly in g/7 
for bounded 7/5. This finishes the proof of the lower bound for Regions 4 and 5. 

If (7/7 — ► as N — > 00, we can use exactly the same strategy, with (4.19) replacing 
the bound (4.18). Considering the case n > 1/x, (4.22) is still valid, but with u(n) 
replaced by 



/ CL\ 1/8 TIT / a \ 

v'(n) = x + const, n y gi + const, n y- J 1 + — y—J 



tit ( a\ V 8 



(4.33) 



For N a > max{2/e, e 2 N} we proceed exactly as above. (We recall that < e < 1, 
and M = eR.) The reason why we have to ensure that N a > e 2 N is the second term 
in (4.33), where we want I to be small. (Note that we choose I = M/(8N a ).) For 
N a < max{2/e, e 2 N}, we replace the bound (4.29) by 

J\f3 at n 

-^efaM/ATJ < -|| m „{2/f, e 2 iV} 

< j^m,)^) ^^, , (4.34, 
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where we used e(x) < and (4.14). Note that, for small g/j, the last fraction is 
order 1. We obtain 

Ne 1 - 



£ QM (iY,L )r , a ) 



with 1Z' given by 

w = 



€ 

1 

X 



X 



e 



> E (N, L, g)1Z' , 



(4.35) 



9 

667 



C 



1-Ce- 



1 + 



e \r/ 

9/1 



i + 



l 



6^(5/7) r 



(4.36) 



2 V~ ' iVeV E LL (l,l,g/-y) 1 + 1 / a 

Here we used again (4.14) to estimate R from above, and e 2 N < N a < N in the 
error term. This proves the lower bound in Region 3, as long as a/(r£(g/j)) stays 
bounded (or at least does not increase too fast as a/r goes to zero with N). Note 
that, for small g/j, 

I a, , u,x,i«.x,i rN 



rt,{gh) r 



£( 7 /^)(«+2)/(«+l) 



with Ltf defined in (2.11), and hence the bound is uniform for rN/L^p bounded. 

As explained briefly in the introduction to this section, if Nr/L^F is not small, 
we have to use Thm. 3.2 instead of Thm. 3.1. Hence we will now assume that 
A = (rjV/Z T F) 1/3 > 1, but still 5/7 -> as N -> 00. Instead of (3.2) we will use the 
bound (3.4) in (4.22). This gives (4.22) with v(ri) replaced by 



v"(n) = C 



7 + x'/» +vv 



+ 



1 + 



Vxr 



na 



(4.38) 



Let < e < 1. For N a > max{2A 2 /e, e 2 iV} we proceed as above, but choosing 
77 = A 2 /(5N a ) for \e < 5 < e. Moreover, we choose x = l/(eN a rf). Using R ~ Ltf, 
this gives (4.35), with 1Z' replaced by 



n" 



l-C 



Na ( e 
iR + [jA 2 



+ 



aA 



ree z 



+ 



1 + 



sl/2 



€€ 



■71/2 




•(l-?) 3 1-Ce 



1 + 



2A*_ 

Ne 2 e 



9h 



(4.39) 



ELL(i,i, 5 / 7 )i+iA 

Note that, for 5/7 < 1, Na/R ~ ap < 1, and (aA/r) 3 = (a/r) 2 Na/R. Moreover, 
,4 2 /iV < 7V-V3 if Z TF /L = (iV £ ,L) 1 /(*+i) i s bounded away from zero. Hence, for 
bounded 1/A, this gives the desired lower bound for Region 3. 

In summary, we have thus shown that, in the limit N — > 00 and r/L — ► 0, 



lim inf 



E® M (N,L,r, a) - Ne^/r 2 



> 1 



(4.40) 



£ LL (/V,L, 5 ) 

uniformly in the parameters, provided (1.12) holds, a/r — > and l/(NgL) stays 
bounded. This finishes the proof of the lower bound for Regions 3-5. 
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4.4 Lower bound for Regions 1 and 2 

We are left with the lower bound for Regions 1 and 2. We proceed similarly to [26, 
Sect. 5.1]. Let pff Lg denote the minimizer of the GP functional (2.4) under the 
normalization condition J K p = N, and let <fi(z) = {p^fi g {z)) 1 / 2 . We write a general 
wave function ^ as 

JV 

*(xi,...,xjv) = F(xi,...,xjv) Y[<f)(z k )b r (-x£) (4.41) 

k=l 

and assume that = 1. In evaluating the expectation value of i?jv,L,r-,a) we use 

partial integration and the GP equation 

-< f) " + V4> + g\4>\ 2 4>=(E GP {N,L 1 g) + 9 - j H 4 )^. (4.42) 

Moreover, we split a fraction of the kinetic energy into a part where the particles are 
closer than a distance T > Rq, and it's complement. This splitting will be important 
in the proof of BEC below. More precisely, for fixed i and Xj, j / i, let 

Xl ' T{ ' \ otherwise , 1 6) 

and let %i,T = 1 — Xi,T- Then, for < e < 1, 

(*\H*,L,r,a\*) = ^r + E GP (N,L,g) + ^j\<P\* + Q(F) 

N 

+(l-e) / |V i F| 2 x i ,r(x i )n^(- z fc) 2ft r(x^) 2 d 3 x fc , (4.44) 
J k=i 

with 

. / N 

= J (E [e|V J F| 2 + (l-e)|V l F| 2 ^,T(x i 

N \ N 

+ ^^(|x,-x,|)|F| 2 - 5 ^|^)| 2 |F| 2 X\<p{z k fb r {*ifd^ k . 

i<j i=l / k=l 

(4.45) 

To bound Q(F) from below, we use again the box method. We divide R in z- 
direction into intervals of length M, labeled by a, put N a particles in box a, neglect 
the interaction between boxes, and minimize over the distribution of the N particles. 
This gives a lower bound. More precisely 

infQ(F) > inf VinfQ Q (F Q ) , (4.46) 

b \ N <xl „ *<* 
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where F a = F Q (xi, . . . , xjv q ), and Q a is the same as Q, but with all the integrations 
restricted to the box a, and N replaced by N a . The infima are under the normaliza- 
tion conditions / \F\ 2 Y\ k= i H^X^if = 1 and J a \F a \ 2 Uk=i <P{z k ) 2 b r {^) 2 = 1, 
respectively. 

We consider two cases. Choose 5 > 0. First, assume that, for all z € a, \4>(z)\ 2 > 
5N/L. We use the same method as in [26, Eqs. (5.28)-(5.34)] to get rid of the <p 2 in 
the measure (f){z) 2 dz. Let </> Q) min and </> Q ,max denote the minimal and maximal value 
of 4> inside the box a, respectively. We first proceed as in (3.30)-(3.36), and obtain, 
for T > r{a/r) 1 ^ = radius of U, 



£ 

i=i 



> 



- e)\ViF\ 2 Xi, T (xi) + \ E ~ x iDI F l 

TO 



\{(p{z k ) 2 b r {^) 2 d^ k 

fc=l 

E / [^^a'C/dxi-x^Dxe^x^/OlFl 2 ] 3 fe) 2 6 r (x^) 2 d 3 x fe . 



Here U is given in (3.33), and a' is given after Eq. (3.36). Denoting 



F(x u . . . ,x Na ) =F(xi,...,Xiv a ) JJ^(zfc) , 



and using 



we get 



|V,F| 2 < 2|V,F| 2 []^,) 2 + 2|F 



k=i 



2 S^Pzec 



k'\2 



k=l 



b 2 ■ 

a, mm 



(4.47) 



(4.48) 



(4.49) 



f \V t F\ 2 f[<p(z k ) 2 b r ( X i) 2 d 3 x k > [ 

Ja k=1 Ja 



Here we denoted 



L 3 

CiV 9 L = — SUp | (/.'(z) | 2 , 



[] 6 r (x^) 2 (i 3 X fc . 
J fc=l 

(4.50) 
(4.51) 



which, by scaling, depends only on NgL. Estimating 4> 2 by it's maximum in the last 
term in (4.45), we therefore have 



Qa(F) > Qa(F) - N a g4>l 



max ' t jj ■> ly ot j 



(4.52) 



with 



Na r ( 6 2 \ N 

Q«( F ) = E / h\^ F \ 2 + 5^-^(1* - **(ol)l*1 2 II b r (*i) 2 d 3 x k . 

i=l • 7a V </>a,max / fc=1 

(4.53) 
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Denote the infimum of Q a (F) over all F (under the normalization condition 
J|F| 2 n fc &r(xfc) 2 (i 3 x A . = 1) by E a (N a ,M), and choose e = (a/r) 1 / 8 . Looking at 
the proof of Thm. 3.1, we see that the lower bound in Subsect. 3.2 was obtained ex- 
actly from an expression like (4.53) (compare with (3.36)), except for the additional 
factor (f)^, min /(pa max- This factor can be estimated by 



9 ■ 



> 1 



L 



2M C NgL 



(4.54) 



Therefore we can apply (3.2), and, in addition, Lemma 2 from Sunsection 3.4 to 
estimate E^P from below. This gives 



L_^^\^-C,(N a ,M)), (4.55) 



where 



v(N,M) = Ny^gM + N 



1/8 



1 + 



Nr /a\i 



M \r 



(4.56) 



Note that v(N, M) is monotone increasing in N. 

This bound is of no use for large N a , however. Therefore we will use again the 
box method, as in Subsect. 4.3 (see also [31]), with small boxes £ = Mr] for some 
rj~ l € N. The use of Lemma 3, with L{x) = x[x — 1] + and A = 4, implies 



1 iV 2 

E a (N a ,M)>- 1 £g[l 



(4.57) 

Let 1 > e > 2/N such that iN £ N, and choose rj = (eiV) -1 . We estimate N a < N 
in the last term in (4.57), and choose M = eL for some e > e. Minimizing over N a 
yields 

E a (N a , M) - g<Pl min N a > -lg<f>i >min MK , (4.58) 

with 



i - c 



egL 1 /ax V8 
e 3 N e \r 



1 



rN /a\V8 



eL \r 



(4.59) 



For boxes a where </> 2 min < NS/L, we just neglect the positive terms in (4.45). 
This gives 

N „ N 



1 



OX V8 

r J 



E QM (N,L,r, a) 
>^ + E^(N,L,g ) + lJ 



E / |V l F| 2 XliT ( Xl )J]0(^) 2 & r (x^) 2 d 3 x fe 
i=i J k=i 



i=i 



+ J] inf £(7V Q , M) - 5 < max iV a 



Q\l/8 JV C NgL 

r) L 2 5 
,N*g 



(4.60) 
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Here we neglected the condition ^ Q N a = N on the N a , which can only lower the 
infimum. Moreover, the error in replacing <^ max by 0^ min in the term in square 
brackets in (4.60) is bounded above by 

M 

2NgM sup 4>{z)W{z)\ = 2N 2 g—C NgL . (4.61) 

Note that, by scaling, C^gL depends only on NgL. 

Using (4.58), J2 a ^a,min M < I 4> A and E a 0a,min M < N and dropping the posi- 
tive second term on the left side of (4.60), we finally obtain 

> -f / |tf(R - I) - Q 1 ' 8 - ffl-i - . (462) 

Note that E GP (N,L,g) > \g f <j) 4 and E GP (N,L,g) > c s N/L 2 for some constant c s 
depending only on s. Hence we see that, for bounded NgL, the parameters e, e and 
5 can be chosen arbitrary small with a/r and AT to show the desired lower bound, as 
long as Nr/L <C (r/a) 1 / 4 and, in particular, for Nr/L bounded. (Note that L ~ L in 
Regions 1 and 2.) Here we also need that C^gL and Cn 9 l are uniformly bounded if 
NgL stays bounded, which follows by the same methods as in the proof of Lemma 1 
in the Appendix. For bigger values of Nr/L we have to proceed differently, using 
the lower bound (3.4) instead of (3.2), as we also did in the previous subsection for 
the lower bound for Region 3. We omit the details. The results of this subsection 
can thus be summarized as 

. c E^ M (N,L,r,a)-Ne ± /r 2 
hm mf E GP (N, L, g) ^ 1 (463) 

in the limit N ^ oo and r/L — > 0, uniformly in the parameters as long as NgL stays 
bounded. This finishes the proof of the lower bound for Regions 1 and 2. 



5 Bose-Einstein Condensation 

In this last section we investigate the question of Bose-Einstein condensation in the 
ground state. It will be proved to occur in Regions 1 and 2 but it probably also 
occurs in part of Region 3; we cannot prove this and it remains an open problem. 
BEC in the ground state means that the one-body density matrix 7(x, x'), which is 
obtained from the ground state wave function by 

7(x,x / ) = N J ^ ( x ,x 2 ,...,x iV )^o(x / ,X2,...,x 7 v)*(i 3 X2---d 3 x 7 v , (5.1) 

factorizes as Nip(x)ip(x.') for some normalized ip (in the N — > oo limit, of course). 
This, in fact, is 100% condensation. It was proved in [25] for a fixed trap potential 
in the Gross-Pitaevskii limit, i.e., for both r/L and Na/L fixed as N — ► oo. Here we 
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extend this result to the case r/L — > with NgL fixed. The function tp is the square- 
root of the minimizer of the ID GP functional (2.4) times the transverse function 

BEC is not expected in Regions 4 and 5. Lenard [21] showed that the largest 
eigenvalue of 7 grows only as N 1 / 2 for a homogeneous gas of ID impenetrable bosons 
and, according to [34] and [10], this holds also for a GT gas in a harmonic trap. (The 
exponent 0.59 in [14] can probably be ascribed to the small number of particles 
(N = 10) considered.) 

Our main result about BEC in the ground state is: 

Theorem 5.1 (BEC in Region 2). If N — ► oo ; r/L — > with NaL/r 2 fixed, then 

±£->y(rx ± ,Lz;rx! ± ,Lz>) - 6(x ± )6(x' ± )^ GP (z)0 GI V) (5.2) 

in trace norm. Here 4> GP is the minimizer of the GP functional (2.4) with N = 1, 
L = 1 and interaction parameter NgL. 

Proof. As in the proof of the energy asymptotics in Region 2, we have to distinguish 
the cases rN/L small or large. For simplicity we consider only the case rN/L <C 
(r/a) 1 / 4 . The case of larger rN/L can be treated in the same manner, replacing the 
bound (3.2) by (3.4) in Subsect. 4.4, and choosing T in (4.43) appropriately. 

From the lower bound to the energy in Subsect. 4.4, together with the upper 
bound in Subsect. 4.1, we infer that if ^0 is the ground state of the Hamiltonian 
HN,L,r,ai T = r(a/r) 1 / 4 and F is defined by 

N 

*o(xi, . . . , xjv) = F(xi, . . . , XJV ) 11 L- 1 / 2 GP (z fe /L)6 r (x^) , (5.3) 



k=i 



then 



I? A " X 



Jim ^ E / I V l F| 2 Xl , T (x,) J] L~ V GP (^/L) 2 6 r (x^) 2 d 3 x fc = (5.4) 
* " 1 ' fc 1 



in the limit N — > 00, r / L — > with NgL fixed. Here XiT is given in (4.43). Note 
that in this limit NT 3 /(r 2 L) = Nr/L(a/r) 3 ^ 4: — > 0, i.e., the volume of the set where 
Xi,T is zero is small compared to the total volume r 2 L. Eq. (5.2) now follows, using 
the methods of [25]. □ 

A Appendix: Proof of auxiliary Lemmas 

Proof of Lemma 1. Without restriction we may assume that V 1 - > 0. The existence, 
uniqueness and positivity of a minimizer <f> p are standard (cf., e.g., [27]). From the 
variational equation 

(_ A ^ + ^(xj.) + 2p|0 p (x ± )| 2 ) M^) = fhM^) (A.l) 
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we infer that, for Kfe 1 - — y -1 ) the integral kernel of (—A ± + 1) 1 , 

M^) = j - y-L) + 1 - ^(y^) - 2p|^(y ± )| 2 ) p (y ± )d 2 y ± • (A.2) 

Using positivity of V ± , 4> p and K, as well as the normalization of <fi p , we obtain the 
bound 

^(x ± )el x± l < sup/ e^^-y^^+lj^/y 1 
x-L Jy ± (y ± )<^ P +l 



< (^p + 1) sup 

x x V- 7v ' ± (y ± )^f+ 1 



1/2 



;l x± lK(x ± -y ± ) 



= C p . 



(A.3) 



Since /_t p is uniformly bounded for p in a bounded interval, so is C p . Moreover, 



|V^ p (x J 



with 



< 



Aip + 1 - F ± (y ± ) - 2p|^(y ± )| 2 ^(y ± )d 2 y ± 

(A.4) 



C p = sup ^ p + 1 - V ± (y ± ) - 2p|0 p (y- L )| :a ^(y^) , 



(A.5) 



which is finite and uniformly bounded for bounded p because of (A.3) and the fact 
that that V 1 - is polynomially bounded at infinity by assumption. □ 



Proof of Lemma 2. We write (1.3) as 

n 

H ™,9 = E~^f + hYl 



■— -d] + g5(z l -z J ) 



(A.6) 



To bound this expression from below, we want to use Temple's inequality, and for 
this purpose we have to smear out the 5-function interaction. We use the following 
Lemma (compare with [2, Lemma 6.3]), whose proof can be found below. 

Lemma 5. Let <9 2 denote the Neumann Laplacian on an interval [0,1], let zq £ 
(0,£), and choose positive numbers A, B and a such that R = B &rctan(BA/2a) < 
mm{zo,£ — zq}. Then 



-ad 2 z + A5(z - zq) - ^9{R - \z - z \) > . 



(A.7) 



We apply this result to the operator in square brackets in (A.6). This gives 
H n , g > ~¥j + \ E ( n -\)B^ {R " lZj ~ ^XW-mM) ' ( A - 8 ) 

.7=1 i^i 
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with R = B arctan(i?<7(n — l)/2) and i? > arbitrary. Temple's inequality (2.61) 
implies that, for 

U( Zl ,...,z n ) = _\\ B 2 6 ( R - \ z i ~ z i\)X[R,£-R])(zi) (A.9) 
and (U k ) = £- n f U k d Zl . . . dz n , 

&M»)*w(i-$ W jhm) • (A ' 10) 

provided the term in the denominator is positive. By Schwarz' inequality, 

(U 2 ) < ^(U) . (A.ll) 

Moreover, 

{u) = T&h Si dz Sr R dw9(R -\ z -^ ) = w^- 2R) - (A - 12) 

Using x > arctan(x) > x(l — x/3) for x > 0, this leads to the estimate 

ciD/ d \ >» i I i\9 f-, B2 9 n B 9 n \ ( n \ 

(A.13) 

under the assumption ir 2 /£ 2 > n(n — l)g/£. Choosing B = i(£gir 2 ) this gives, 
for K = {n 2 lgf/ 2 /^, 

<(n, i, g) > \n(n - \)\ (l - U - ^^ 2 ) (l - ^) (A.14) 

and proves the desired lower bound. □ 
Proof of Lemma 5. Let h denote the operator on the left side of (A. 7). The function 

f( z ) = S cos \z-z \)/B^j for < \z - z \ < R 

^ 1 otherwise 

fulfills the Schrodinger equation hf = 0, and since it is positive, it must be the 
ground state of h. Hence h > 0. □ 

Proof of Lemma 3. Let p = (XN/M). For n < p we use the estimate 

E(n) > L{n)K(p) , (A.16) 

which follows from monotonicity of K. For n > p, we use superadditivity, which 
implies that 

Tl 

E(n) > [n/p]E(p) > —L(p)K(p) , (A.17) 
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where [x] denotes the integer part of x. Now let t = J2n<p c « n - By convexity of L 
and the fact that L(0) = 0, 

^c n L(n) > ML(t/M) . (A.18) 

n<p 



Hence 

n>0 



ML(t/M) + (N -t)- 



2p 



(A.19) 



To obtain a lower bound, we have to minimize the right side over all < t < N. 
Note that (A.19) is convex in t, hence the minimum is either taken at t = N, or, if 
it is taken at some to < N, it's right derivative at to has to be positive. Using (3.83) 
this leads to 

Ij&m > M . (A.20) 
2Xt /M ~ 2p y J 

By our assumptions on L, L(x)/x is monotone increasing, and hence to > pM/X > N. 
Thus we can set t = N in (A.19), and obtain the desired lower bound. □ 
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